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^ ■ Abstract 

An exploratory study of charmless 3-body decays of B mesons is presented using a simple model 
based on the framework of the factorization approach. The nonresonant contributions arising from 
B —* -P1-P2 transitions are evaluated using heavy meson chiral perturbation theory (HMChPT). 
The momentum dependence of nonresonant amplitudes is assumed to be in the exponential form 

■ e ~ a ]siRPB-{pi+Pj) so that the HMChPT results are recovered in the soft meson limit pi, pj — > 0. In 
addition, we have identified another large source of the nonresonant signal in the matrix elements 

t-h ■ of scalar densities, e.g. (KK\ss\0), which can be constrained from the decay B — > KsKsKs 

■"sj- ■ or B~ —> K~ K$Ks- The intermediate vector meson contributions to 3-body decays are identified 

through the vector current, while the scalar meson resonances are mainly associated with the scalar 
density. Their effects are described in terms of the Breit-Wigner formalism. Our main results are: 

(i) All KKK modes are dominated by the nonresonant background. The predicted branching 
^ . ratios of K + K~K S ( L j, K + K~ K~ and K~KsKs modes are consistent with the data within errors. 

(ii) Although the penguin-dominated B° — » K + K~K$ decay is subject to a potentially significant 
tree pollution, its effective sin 2(5 is very similar to that of the KsK$Ks mode. However, direct 
CP asymmetry of the former, being of order —4%, is more prominent than the latter, (iii) For 
B — > Kirn decays, we found sizable nonresonant contributions in K~Tr + ir~ and ~K°i: + -n~ modes, 
in agreement with the Belle measurements but larger than the BaBar result, (iv) Time-dependent 
CP asymmetries in Kg^ °vr , a purely CP-even state, and Ks^ 7r , an admixture of CP-even 
and CP-odd components, are studied, (v) The 7r + 7r~7r° mode is found to have a rate larger than 
7r + 7r~7r~ even though the former involves a ir° in the final state. They are both dominated by 
resonant p contributions, (vi) We have computed the resonant contributions to 3-body decays 
and determined the rates for the quasi-two-body decays B — > VP and B — > SP. The predicted 
pn, fo(980)K and /o(980)7r rates are in agreement with the data, while the calculated (pK, K*n, pK 
and Kq (1430)-7r are in general too small compared to experiment, (vii) Sizable direct CP asymmetry 
is found in K + K~ K~ and K + K~tt~ modes. 
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I. INTRODUCTION 



Three-body decays of heavy mesons are more complicated than the two-body they 
receive resonant and nonresonant contributions and involve 3-body matrix elements. The three- 
body meson decays are generally dominated by intermediate vector and scalar resonances, namely, 
they proceed via quasi-two-body decays containing a resonance state and a pseudoscalar meson. 
The analysis of these decays using the Dalitz plot technique enables one to study the properties 
of various resonances. The nonresonant background is usually believed to be a small fraction of 
the total 3-body decay rate. Experimentally, it is hard to measure the direct 3-body decays as the 
interference between nonresonant and quasi-two-body amplitudes makes it difficult to disentangle 
these two distinct contributions and extract the nonresonant one. 

The Dalitz plot analysis of 3-body B decays provides a nice methodology for extracting infor- 
mation on the unitarity triangle in the standard model (SM). For example, the Dalitz analysis 
combined with isospin symmetry allows one to extract the angle a from B — > pit — > -nitir 
Recently, a method has been proposed in [2| for determining CKM parameters in 3-body decays 
B — > Ktttt and B s — » K-kit. This method was extended further in [3] to AI = l,I(K*it) = 1/2 
amplitudes in the above decays and to / = 1 amplitudes in B s — > K*K and B s — » K*K. 

Nonresonant 3-body decays of charmed mesons have been measured in several channels and the 
nonresonant signal in charm decays are found to be less than 10% [4]. In the past few years, some 
of the charmless B to 3-body decay modes have also been measured at B factories and studied 
using the Dalitz plot analysis. The measured fractions and the corresponding branching ratios of 
nonresonant components for some of 3-body B decay modes are listed in Table HI We see that 
the nonresonant 3-body decays could play an essential role in B decays. It is now well established 
that the B — ► KKK modes are dominated by the nonresonant background. For example, the 
nonresonant fraction is about 90% in W -> K+K-K decay. While this is a surprise in view of 
the rather small nonresonant contributions in 3-body charm decays, it is not entirely unexpected 
because the energy release scale in weak B decays is of order 5 GeV, whereas the major resonances 
lie in the energy region of 0.77 to 1.6 GeV. Consequently, it is likely that 3-body B decays may 
receive sizable nonresonant contributions. At any rate, it is important to understand and identify 
the underlying mechanism for nonresonant decays. 

The direct nonresonant three-body decays of mesons in general receive two distinct contributions: 
one from the point-like weak transition and the other from the pole diagrams that involve three- 
point or four-point strong vertices. For D decays, attempts of applying the effective SU(4) x SU(4) 
chiral Lagrangian to describe the DP — ► DP and PP — > PP scattering at energies ~ mo have 
been made by several authors [l^i (3) 20, 21, 22] to calculate the nonresonant D decays, though 



21 



221, the 



in principle it is not justified to employ the SU(4) chiral symmetry. As shown in 
predictions of the nonresonant decay rates in chiral perturbation theory are in general too small 
when compared with experiment. With the advent of heavy quark symmetry and its combination 



with chiral symmetry 23|, |24j, |25j , the nonresonant D decays can be studied reliably at least in the 



kinematic region where the final pseuodscalar mesons are soft. Some of the direct 3-body D decays 



were studied based on this approach [26j, |27|] . 
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TABLE I: Branching ratios of various charmless three-body decays of B mesons. The fractions 
and the corresponding branching ratios of nonresonant (NR) components are included whenever 



available. The first and second entries are BaBar and Belle results, respectively. 

Decay BR(1(T 6 ) BR NR (l(r 6 ) NR fraction(%) Ref. 

B-^tt+tt-tt- 16.2 ±1.2 ±0.9 2.3 ± 0.9 ± 0.5 < 4.6 13.6 ±6.1 [5j 

B K-TT+TT- 64.1 ±2.4 ±4.0 2.87 ± 0.65 ± 0.43±g;|| 4.5 ±1.5 [6j 

48.8 ±1.1 ±3.6 16.9 ± 1.3 ± 1.3lJ;5 34.0 ± 2.2^-i [7J 

B-^K + K-K~ 35.2 ± 0.9 ± 1.6 a 50 ± 6 ± 4 141 ± 16 ± 9 [8] 

32.1 ± 1.3 ±2.4 b 24.0 ± 1.5 ± 1.5 c 74.8 ± 3.6 c [9j 

B-^K-K S K S 10.7 ±1.2 ±1.0 [10] 

13.4 ±1.9 ±1.5 [11] 
W K u ir + ir- 43.0 ±2.3 ±2.3 [12]" 

47.5 ±2.4 ±3.7 19.9 ± 2.5 ± 1.6±?J 41.9 ± 5.1 ± 0.6±|;| [13] 
B° -> K-tt+tt 34.9 ±2.1 ±3.9 < 4.6 [14] 

36.6ita ± 3.0 S^lf^f < 9.4 [15] 

Tf^K+K-K 23.8 ±2.0 ± 1.6 26.7 ± 4.6 112.0 ± 14.9 [16] 

28.3 ±3.3 ±4.0 [11] 

5° - Er 5 ^5^5 6.9l°| ± 0.6 [17] 

4.2±};|±0.8 [11] 



"When the intrinsic charm contribution is excluded, the charmless branching ratio will become (33.5 ± 0.9 ± 
1.6) x 1CT 6 . 

b When the contribution from B + — > \cqK + is excluded, the charmless branching ratio will become (30.6 ± 
1.2 ±2.3) x 10~ 6 . 

c Belle found two solutions for the fractions and branching ratios. We follow Belle to use the large solution. 



For the case of B mesons, consider the three-body B decay B — > P1P2P3. Under the factorization 
hypothesis, one of the nonresonant contributions arises from the transitions B — > P\Pi- The 
nonresonant background in charmless three-body B decays due to the transition B — > P1P2 has 
been studied extensively 



28, 



29, 



30, 



31 



32, 



331 ] based on heavy meson chiral perturbation theory 



(HMChPT) [23j, [24j, [25|] . However, the predicted decay rates are, in general, unexpectedly large. For 
example, the branching ratio of the nonresonant decay B~ — > 7r + 7r~7r - is predicted to be of order 
10~ 5 in 28|, 129J, which is too large compared to the limit 4.6 x 10~ 6 set by BaBar [5j. Therefore, 



it is important to reexamine and clarify the existing calculations. 

The issue has to do with the applicability of HMChPT. In order to apply this approach, two 
of the final-state pseudoscalars in B — > P1P2 transition have to be soft. The momentum of the 
soft pseudoscalar should be smaller than the chiral symmetry breaking scale of order 1 GeV. For 
3-body charmless B decays, the available phase space where chiral perturbation theory is applicable 
is only a small fraction of the whole Dalitz plot. Therefore, it is not justified to apply chiral and 
heavy quark symmetries to a certain kinematic region and then generalize it to the region beyond 
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its validity. In this work we shall assume the momentum dependence of nonresonant amplitudes in 
the exponential form e~ a ^ PB '^ Pi+p ^ so that the HMChPT results are recovered in the soft meson 
limit pi, pj — > 0. We shall see that the parameter q nr can be fixed from the tree-dominated decay 
B~ — > ir + 7r~ir~ . 

However, the nonresonant background in B — > P1P2 transition does not suffice to account for 
the experimental observation that the penguin-dominated decay B — ► KKK is dominated by the 
nonresonant contributions. This implies that the two-body matrix element e.g. (KK\ss\0) induced 
by the scalar density should have a large nonresonant component. In the absence of first-principles 
calculation, we will use the 1?° — ► KgKgKg mode in conjunction with the mass spectrum in 
B° — > K + K~K° to fix the nonresonant contribution to (KK\ss\0). 

In this work, we shall study the charmless 3-body decays of B mesons using the factorization 
approach. Besides the nonresonant background as discussed above, we will also study resonant 
contributions to 3-body decays. Vector meson and scalar resonances contribute to the two-body 
matrix elements (-Pii^l^lO) and (PiP2\S\0), respectively. They can also contribute to the three- 
body matrix element (-Pii^lKt — A^B). Resonant effects are described in terms of the usual 
Breit-Wigner formalism. In this manner we are able to figure out the relevant resonances which 
contribute to the 3-body decays of interest and compute the rates of B — » VP and B — > SP. 
In conjunction with the nonresonant contribution, we are ready to calculate the total rates for 
three-body decays. 

It should be stressed from the outset that in this work we take the factorization approximation 
as a working hypothesis rather than a first-principles starting point. If we start with theories such 
as QCD factorization 34 |. or pQCD [35(] or soft-collinear effective theory 36f] , then we can take 
power corrections seriously and make an estimation. Since factorization has not been proved for 
three-body B decays, we shall work in the phenomenological factorization model rather than in the 
established theories such as QCDF. That is, we start with the simple idea of factorization and see 
if it works for three-body decays, in the hope that it will provide a useful zeroth step for others to 
try to improve. 

The penguin-induced three-body decays B° — ► K + K~Ks and KgKgKg deserve special atten- 
tion as the current measurements of the deviation of sin2/3 e fj in KKK modes from sin2/3j/^ g 
may indicate New Physics in b — ► s penguin-induced modes. It is of great importance to examine 
and estimate how much of the deviation of sin2/3 e fj is allowed in the SM. Owing to the presence 
of color-allowed tree contributions in B° — > K + K~ K$, this mode is subject to a potentially signif- 
icant tree pollution and the deviation of the mixing-induced CP asymmetry from that measured 
in B — > J/ipKs could be as large as O(0. 10). Since the tree amplitude is tied to the nonresonant 
background, it is very important to understand the nonresonant contributions in order to have a 
reliable estimate of sin2/3 e g- in KKK modes. 

The layout of the present paper is as follows. In Sec. II we shall apply the factorization 
approach to study B° — ► K + K~K$ and KsK$Ks decays and discuss resonant and nonresonant 
contributions. In order to set up the framework for calculations we will discuss B — » KKK modes 
in most details. We then turn to Ktttt modes in Sec. III. The tree-dominated modes KKtt in 
Sec. IV, and tttttt in Sec. V. In Sec. VI, we determine the rates for B — > VP and B — ► SP and 
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compare our results with the approach of QCD factorization. Sec. VII contains our conclusions. 
The factorizable amplitudes of various B — > P1P2P3 decays are summarized in Appendix A. The 
relevant input parameters such as decay constants, form factors, etc. are collected in Appendix B. 



II. B -> KKK DECAYS 



For 3-body B decays, the b — ► sqq penguin transitions contribute to the final states with odd 
number of kaons, namely, KKK and Kirir, while b — > uqq tree and b — * dqq penguin transitions 
contribute to final states with even number of kaons, e.g. KKir and tttttt. We shall first discuss the 
b — * s penguin dominated 3-body decays in details and then turn to b — > u tree dominated modes. 
For B — > KKK modes, we shall first consider the neutral B decays as they involve mixing- induced 
CP asymmetries. 



A. B° — ► KKK decays 

We consider the decay B° — ► K + K~K° as an illustration. Under the factorization approach, 
the B° — > K + K~K° decay amplitude consists of three distinct factorizable terms: (i) the current- 
induced process with a meson emission, (B° — > K + K°) x (0 — > K~), (ii) the transition process, 
(B° -> K°) x (0 -» K + K~), and (iii) the annihilation process (B° -» 0) x (0 -» K + K~T<°), where 
(.A — > B) denotes ai->B transition matrix element. In the factorization approach, the matrix 
element of the B — > K KK decay amplitude is given by 



(KKK\H eS \B) = % E >i a) (KKK\T p \B), (2.1) 



(s) 

where Ap = V p bV* s and the explicit expression of T p in terms of four-quark operators is given in 
Eq. (|A2p . The factorizable i?° — ► decay amplitude is given in Eq. (|A4|) . Note that 

the OZI suppressed matrix element {K + K~\(dd)v-A\0) is included in the factorizable amplitude 
since it could be enhanced through the long-distance pole contributions via the intermediate vector 
mesons such as p° and uj. Likewise, the OZI-suppressed matrix elements (K + K~\(db) v _ A \B°) and 
(K + K~\d(l — ^y^)b\B°) are included as they receive contributions from the scalar resonances like 
/ (980). _ 

For the current-induced process, the two-meson transition matrix element (K°K Jr \(ub)v-A\B ) 



has the general expression 37] 



(K (pi)K + (p 2 )\(ub)v~A\B ) = ir(p B - Pi + i^+(P2 +Pl)n + w-(P2 ~ Pi)n 

+he pva pp u B (p2+PiT{p2-piY, (2.2) 

where {qiq 2 ) v _ A = fiT/LtC 1 _ 75)^2- This leads to 

^CSL = (K-(p 3 )\(s U ) v ^A\0)(K\ Pl )K + (p2)\(ub) v ^A\B ) 

= 2m\r + (m| - s i2 - m\)u+ + (s 2 3 - S13 - m| + mi)w_ , (2.3) 
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if k" 



K- 



K- 



(d) 



FIG. 1: Point-like and pole diagrams responsible for the 23° — > K + K v matrix element induced by 
the current ^7^(1 — 75)^, where the symbol • denotes an insertion of the current. 



where Sy = (pj + Pj) 2 ■ To compute the form factors r, u>± and h, one needs to consider not 
only the point-like contact diagram, Fig. 1(a), but also various pole diagrams depicted in Fig. 1. 
In principle, one can apply HMChPT to evaluate the form factors r, uj + and uj- 37j]. However, 
this will lead to too large decay rates in disagreement with experiment [38j. The heavy meson 
chiral Lagrangian given in 



23 



24 



25|] is needed to compute the strong B*BP, B*B*P and BBPP 



vertices. The results for the form factors are 
g fB* s m B ,^/m B mBi 



29 



37] 



UJ+ 



K s 23 - m Ba 
g fB*rnB*y/m B mB 
fl s 23 
Ib _ [b_ 



m 



B* 

PB ■ (j>2 



fn (PB -P1-P2 

m B m B * 

ml 




Pi) - Pi 



S23 



m 



B* 



Pi)P2-{PB -Pi)/m%, 



(2.4) 



fw {PB - Pi- P2) 2 ~ m 2 B s 2 3 - m B « 

where f n = 132 MeV, g is a heavy-flavor independent strong coupling which can be extracted from 
the CLEO measurement of the D* + decay width, \g\ = 0.59 ± 0.01 ± 0.07 [39]. We shall follow [23(] 
to fix its sign to be negative. The point-like diagram Fig. 1(a) characterized by the term f B /(2f%) 
contributes to the form factors u+ and r, while Figs. 1(b) and 1(d) contribute to r and Fig. 1(c) 
contributes to all the form factors. 

A direct calculation indicates that the branching ratio of B° —> K + K~K° arising from the 
current-induced process alone is already at the level of 77 x 10~ 6 which exceeds the measured 
total branching ratio of 25 x 10 -6 (see Table |T]). The issue has to do with the applicability of 
HMChPT. In order to apply this approach, two of the final-state pseudoscalars (K + and K° in 
this example) have to be soft. The momentum of the soft pseudoscalar should be smaller than 
the chiral symmetry breaking scale A x of order 0.83 — 1.0 GeV. For 3-body charmless B decays, 
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the available phase space where chiral perturbation theory is applicable is only a small fraction of 
the whole Dalitz plot. Therefore, it is not justified to apply chiral and heavy quark symmetries 
to a certain kinematic region and then generalize it to the region beyond its validity. If the soft 
meson result is assumed to be the same in the whole Dalitz plot, the decay rate will be greatly 
overestimated. 

In 38|, |40|] we have tried to circumvent the aforementioned problem by applying HMChPT only 



to the strong vertex and use the form factors to describe the weak vertex. Moreover, we introduced 
a form factor to take care of the off-shell effect. For example, Fig. 1(c) can be evaluated by 
considering the strong interaction B~° — * K°B* S followed by the weak transition B* — > K + and the 
result is 

A 



l Fig.l(c) 



Ik 9y/m B m B i 



F(s 2 3,m B *)F*° K (m 



h s 23 



m 



LSI 



VtlB + 



■S23 

m B 



rn 



m B - 



B 



•S23 



F*° K {mi) 



2 i 2 

rrii + m 3 



Sl3 + 



(s 2 3 -m,2 + ra 2 )(m| - s 2 3 - m\ ) 



2m 2 



(2.5) 



where Fq{ k are the B s 



K weak transition from factors in the standard convention 



411 ] and we 



have introduced a form factor F(s23,m B *) to take into account the off-shell effect of the B* pole [401 ] . 
It is parameterized as F(s23,m B *) = (A 2 



m 



LS. 



) /(A 2 — S23) with the cut-off parameter A chosen to 
be A = m B * + Aqcd • Needless to say, this parametrization of the form factor is somewhat arbitrary. 
Moreover, the nonresonant contribution thus calculated is too small compared to experiment. 

The Dalitz plot analysis of — > K + K^T^ has been recently performed by BaBar 16]. In the 
BaBar analysis, a phenomenological parametrization of the non-resonant amplitudes is described 
by 



A 



NR 



)(1 + ^NRe 



i(/3+<5 NR )^ 



and resonant terms are described by 
A R = Y / Cr(l + b r )f r e^ 



-+<5r+/3) 



A R = ^c r (l - b r )f r e 



i(</>r-S r +f3) 



(2.6) 



(2.7) 



The BaBar results for isobar amplitudes, phases and fractions from the fit to the B° — ► K + K~K° 
are summarized in Table [TIJ It is evident that this decay is dominated by the nonresonant back- 
ground. For our purpose, we will parametrize the current-induced nonresonant amplitude Eq. (|2,3p 

as 

A HMChPT 



A r 



;PS-(P1+P2) e i</>12 



(2.8 



i current— ind ^-current— ind e 

so that the HMChPT results are recovered in the chiral limit pi, P2 — > 0. That is, the nonresonant 
amplitude in the soft meson region is described by HMChPT, but its energy dependence beyond 
the chiral limit is governed by the exponential term e~ Q NR PB '( pi +P 2 ). In what follows, we shall use 
the tree-dominated B~ — ► 7r + 7r~7r~ decay data to fix a NR , which turns out to be 



0.103; 



-0.018 
-0.01 



*GeV- 



(2.9) 



^ 0(l/(2m B A x )) based on the dimensional 
argument. The phase 4>i2 of the nonresonant amplitude in the {K + K ) system will be set to zero 



This is very close to the naive expectation of 
argument. Th 
for simplicity. 
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TABLE II: BaBar results for isobar amplitudes, phases, and fractions from the fit to the B° — ► 
K + K~K° Three rows for non-resonant contribution correspond to coefficients of exponential 
functions in Eq. (|2.6p . while the fraction is given for the combined amplitude. For the nonresonant 
decay mode in K + K~, the amplitude c\2 and the phase <^>i2 in Eq. (|2.6|) are fixed to be one and 
zero, respectively. Errors are statistical only. 



Decay 


Amplitude c r 


Phase 4> r 


Fraction (%) 


0(1O2O)K° 


0.0085 ± 0.0010 


-0.016 ±0.234 


12.5 ± 1.3 


Zo(980)if° 


0.622 ± 0.046 


-0.14 ±0.14 


40.2 ± 9.6 


X (1550)X° 


0.114 ±0.018 


-0.47 ±0.20 


4.1 ± 1.3 


(K+K~) NR K° 


1 (fixed) 


(fixed) 




(K+K°) NR K- 


0.33 ± 0.07 


1.95 ±0.27 


112.0 ± 14.9 


(K-K°) NR K+ 


0.31 ±0.08 


-1.34 ±0.37 




Xco(lP)K° 


0.0306 ± 0.00649 


-I'M ± 0-54 


3.0 ± 1.2 


D+K" 


1.11 ±0.17 




3.6 ±1.5 


D+K- 


0.76 ±0.14 




1.8 ±0.6 



For the transition amplitude, we need to evaluate the 2-kaon creation matrix element which can 
be expressed in terms of time-like kaon current form factors as 

7K+K- 



(K + (p K+ )K-(p K -)\q 7fl q\0) 



) F 



(PK+ ~ Pk- 

(K° (p K0 )K° (p R0 )\ q ltl q\0) = (p K o -p R0 )pt ^ ' ' h \ (2.10) 

The weak vector form factors F^ +K and F^°^-° can be related to the kaon electromagnetic (e.m.) 
form factors F^ K and F^^° for the charged and neutral kaons, respectively. Phenomenologi- 
cally, the e.m. form factors receive resonant and nonresonant contributions and can be expressed 
by 



? K+K~ 



F p + F u + F ()} + F NR , 



? K°K° 



-F + F LU + Fs + Fi 



NR- 



(2.11) 



It follows from Eqs. OTTUD and OTP) that 



K+K- 



K+K- 



F. 



K°K° 
u 

K°K° 



F p + 3F U + -(3F NR 
—F p ± 3-F^, 

-3^-^(3^ + 2^ 



NR)' 



(2.12) 



where use of isospin symmetry has been made. 

The resonant and nonresonant terms in Eq. (|2. 1 1 j) can be parametrized as 



F h (s 23 ) 



Ch 



s 2 3 - im h r h ' 



NR 



S23 



'23 



ln|^ 

A 2 



(2.13) 



with A ~ 0.3 GeV. The expression for the nonresonant form factor is motivated by the asymptotic 
constraint from pQCD, namely, F(t) — ► (l/i)[ln(t/A 2 )] -1 in the large t limit 42|. The unknown 
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parameters Ch,Xi and x\ are fitted from the kaon e.m. data, giving the best fit values (in units of 
GeV 2 for c h ) 0: 



(2.14) 



and 



c p = 3c w = c^ = 0.363, c p(1 4 50 ) = 7.98 x 10 3 , c p( i 700 ) = 1.71 x 10 3 , 
=o;(i420) = -7-64 x 10~ 2 , c w(1650 ) = -0.116, c 0(168O ) = -2.0 x 10~ 2 , 

xi = -3.26 GeV 2 , x 2 = 5.02 GeV 4 , x' x = 0.47 GeV 2 , x' 2 = 0. (2.15) 



Note that the form factors F pw a in Eqs. (|2.11|) and (|2.12|) include the contributions from the vector 
mesons p(770), p(1450), p(1700), w(782), w(1420), w(1650), 0(1020) and 0(1680). It is interesting 
to note that (i) the fitted values of cy are very close to the vector meson dominance expression 
9v 9vkk for V = p, uj,4> 0, 44], where g Vi is the e.m. coupling of the vector meson defined by 
(^|iem|0) = gv^v and 9VKK is the V -> KK strong coupling with -g^K+K- - 9 p k+k-/V^ = 
9luK+k- A/2 — 3.03, and (ii) the vector-meson pole contributions alone yield F^ K (0) m 1,-1 
and F^ +K (0) ~ as the charged kaon does not contain the valence d quark. The matrix element 
for the current-induced decay process then has the expression 

(K ( Pl )\(sb) v ^ A \W)(K + (p 2 )K-(p 3 )\(q q ) v . A \0) = (s 12 - s 13 )Ff K (s 23 )F q K+K ~ (s 23 ). (2.16) 

We also need to specify the 2-body matrix elements (K + K~\ss\0)(K°\sb\B ) induced from the 
scalar densities. The use of the equation of motion leads to 

(K°(p 1 )\sb\B (p B )) = ^-< F( f*( S23) , (2.17) 

m b - m s 

The matrix element {K + K~\ss\0} receives resonant and non-resonant contributions: 

(k+( P2 )k-(p 3 )\ss\o) =. fr x ~M=E 7 0ifkgf0i ^ K+K r +f?*, 

fs R = 1(3Fnr + 2F' NR ) + - im e-°"», (2.18) 

where /oj denote the generic /o-type scalar mesons, /oj = /o(980), /o(1370), /o (1500), Xq (1550), • • •, 
the scalar decay constant /| Q is defined by (/oj|ss|0) = mf 0i ff . [see Eq. (|Bip ]. gf°i~* K K is the 
/oj — ► K + K~ strong coupling, and the nonresonant terms are related to those in F^ +K through 
the equation of motion. The presence of the nonresonant <7 NR term will be explained shortly. The 
main scalar meson pole contributions are those that have dominant ss content and large coupling 
to KK. We consider the scalar mesons /o(980) and Xo(1550) (denoted as /x(1500) by Belle) which 
are supposed to have the largest couplings with the KK pair. Note that the nature of the broad 
state Xo(1550) observed by BaBar and Belle, for example, what is its relation with /o(1500), is 
not clear. To proceed with the numerical calculations, we shall use gfo( 9SO )^ K+ K = 4.3 GeV, 1 



1 This is different from the coupling gf°( 980 )-^ K+K = 1.5 GeV originally employed in [40j. The coupling 
g/o(980)-+7r+7r ^ -^ 23 GeV can be fixed from a recent Belle measurement of T(/o(980) — ► tt + tt~) [see Eq. 
([535]) ]. Using the BES result ( g MaBO)^KK / g Mf>no)^™)2 = 421 ± 0.25 ± 0.21 Q, one can deduce that 
g f (980)^KK = 2.7 ± 0.6 GeV. In this work, we found that a slightly large coupling g h{mo)^KK wffl giye 
better numerical results. 
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gX (i5BO)^K+K- = 1A GeV , r /o(980) = 80 MeV, r Xo(1550) = 0.257 GeV fl, f Mm (v = m b /2) ~ 
0.46 GeV 46] and // (i53o) — 0.30 GeV. The sign of the resonant terms is fixed by f^ +K (0) = v 
from a chiral perturbation theory calculation (see, for example, [47]). It should be stressed that 
although the nonresonant contributions to and F^ K are related through the equation of 

motion, the resonant ones are different and not related a priori. As stressed in 40(], to apply the 
equation of motion, the form factors should be away from the resonant region. In the presence 
of the resonances, we thus need to introduce a nonresonant term characterized by the parameter 
ctnr in Eq. f|2. 18j) which will be specified later. The parameter a appearing in the same equation 
should be close to the value of a NR given in Eq. (|2.9p . We will use the experimental measurement 
a = (0.14 ± 0.02) GeV" 2 Q. 

As noticed before, the matrix elements (K + K~\(db) v _ A \B ) and (K + K~\d(l — ^)b\B 
included in Eq. (|A4p as they receive intermediate scalar pole contributions. More explicitly, 

J ^K+K- 



arc 



(K + (p 2 )K-(p 3 )\(db) v _ A \BY = El2^ — — (falWy-Jg" 



m 



fOi 



S23 ~ im fo r fo 



(2.19) 



Hence, 



(K v (p 1 )\(sd) v _ A \0)(K + ( P2 )K~(p 3 )\(db\ 



\B°) R 



E 



g J 



f K F^(m 2 K )(m 



B 



2 > 



(2.20) 



The superscript u of the form factor F ' reminds us that it is the uu quark content that gets 
involved in the B to fo i form factor transition. In short, the relevant /o(980) pole contributions to 



B° -> K + K~K° are 



(K°K + K-\T P \B°) 



fa 
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f (980)^K+K- 



m 



im fo T fo 



, m fo 



f%F BK (m%)(m 



m 2 K ) 



2 a ' 



fKF^im^imj 



1 v 
2 ai0 



-a 
2 



(2.21) 



where we have employed Eq. (|2.18p and applied equations of motion to the matrix elements 
(K°\sj5d\0){K + K^\d^/^b\B°) . Comparing this equation with Eq. (A6) of 48], we see that the 
expression inside {• • •} is identical to that of 5° — > /o(980)K°, as it should be. 

We digress for a moment to discuss the wave function of the /o(980). What is the quark structure 
of the light scalar mesons below or near 1 GeV has been quite controversial. In this work we shall 
consider the conventional qq assignment for the /o(980). In the naive quark model, the flavor wave 
functions of the /o(980) and <r(600) read 

1 , 



a 



V2 



uu - 



dd), 



fo = ss, 



(2.22) 



where the ideal mixing for /o and a has been assumed. In this picture, /o(980) is purely an ss 
state. However, there also exist some experimental evidences indicating that /o(980) is not purely 
an ss state. First, the observation of T(J/tp — > fou) « ^T(J/tp — > /o^) 0] clearly indicates the 
existence of the non-strange and strange quark content in /o(980). Second, the fact that /o(980) 
and ao(980) have similar widths and that the /o width is dominated by irir also suggests the 
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composition of uu and dd pairs in /o(980); that is, /o(980) — > tttt should not be OZI suppressed 
relative to ao(980) —* irr]. Therefore, isoscalars o"(600) and /o must have a mixing 

|/ (980)) = | ss) cos6 + | nn) sinfl, k(600)) = -\ss) sin 9 + \nn) cos 9, (2.23) 

with nn = (uu + dd) / \/2. Experimental implications for the /q— a mixing angle have been discussed 



in detail in [49J]. It is found that 9 lies in the ranges of 25° < 9 < 40° and -40° < 9 < -15° (or 
140° < 9 < 165°). Note that the phenomenological analysis of the radiative decays (j) —* /o (980)7 
and _/o(980) — > 77 favors a solution of the 9 to be negative (or in the second quadrant). In this 
work, we shall use 9 = —25°. 

Finally, the matrix elements involving 3-kaon creation are given by 38f] 

(K ( Pl )K + (p2)K-(p 3 )\(M) v ^ A \0)(0\(db) v ^ A \B ) « 0, (2.24) 

(^ (pi)^ + fe)^-(P3)N-75d|0)(0|J756|B ) = v S -f&- { 1 - fi i3-^-m§ \ F X^ (m | ); 



/„.m 6 \ m 2 B -m 



K 



where 



(2.25) 

m u + m s m s -m d ' 

characterizes the quark-order parameter (qg) which spontaneously breaks the chiral symmetry. 
Both relations in Eq. (12.241) are originally derived in the chiral limit 38J] and hence the quark 
masses appearing in Eq. (|2.25|) are referred to the scale ~ 1 GeV . The first relation reflects 
helicity suppression which is expected to be even more effective for energetic kaons. For the second 
relation, we introduce the form factor F KKK to extrapolate the chiral result to the physical region. 
Following H we shall take F KKK (q 2 ) = 1/[1 - (q 2 /A 2 x )} with A x = 0.83 GeV being a chiral 
symmetry breaking scale. 

To proceed with the numerical calculations, we need to specify the input parameters. The 
relevant CKM matrix elements, decay constants, form factors, the effective Wilson coefficients 
a? and the running quark masses are collected in Appendix B. As for the parameter cr NR in Eq. 
(|2.18p . in principle we can set its phase 4>n to zero and use the measured KsK$Ks rate, namely, 



B(B -> K S K S K S ) = (6.2±0.9) x 10 -6 [50|, to fix the parameter ct nr and then use the data obtained 



from the Dalitz plot analysis to determine the strong phases 4> r for resonant amplitudes. However, 
in doing so one needs the data of invariant mass spectra. In the absence of such information, instead 
we will treat (j) a as a free parameter and do not assign any other strong phases to the resonant 
amplitudes except for those arising from the Breit-Wigner formalism. It turns out that if (j) a is 
small, the K + K~ mass spectrum in B° — > K + K~ K$ will have a prominent hump at the invariant 
mass rrix+K- = 3 GeV, which is not seen experimentally (see Fig. 12(c)). We found that (j) a pa 7r/4 
will yield K + K~ mass spectrum consistent with the data 

^ NR = e w/4 (3.36±£$) GeV. (2.26) 



Note that the phase of u NR is consistent with the BaBar measurement shown in Table HH namely, 

^BaBar = L l 9 ± Q.37. 

The calculated branching ratios of resonant and nonresonant contributions to B° — ► K + K~K° 
are summarized in Table IIIII The theoretical errors shown there are from the uncertainties in (i) 
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TABLE III: Branching ratios (in units of 10 -6 ) of resonant and nonresonant (NR) contributions to 
B° — » K + K~ K° . Theoretical errors correspond to the uncertainties in (i) a NR , (ii) m s , F^ K and 
o" NR , and (hi) 7 = (59 ± 7)°. We do not have l/m& power corrections within this model. However, 
systematic errors due to model dependent assumptions may be sizable and are not included in the 
error estimates that we give. Experimental results are taken from Table ILT1 

Decay mode BaBar [16] Theory 

W 2.98 ± 0.45 2.6±g : 8±8:t±8:8 

/o(980)K° 9.57 ±2.51 5.8±S;g±gi±S:8 

X Q {lbW)K~ 0.98 ± 0.33 0.93±8;88±g;ii±8:88 

NR 26.7 ±4.6 18.1±8;|±|;|±g;| 

total 23.8 ± 2.0 ± 1.6 IQ-Stattojtoi 



the parameter a^R which governs the momentum dependence of the nonresonant amplitude, (ii) 
the strange quark mass m s , the form factor Fq K and the nonresonant parameter <r NR , and (hi) the 
unitarity angle 7. 

In QCD calculations based on a heavy quark expansion, one faces uncertainties arising from 
power corrections such as annihilation and hard-scattering contributions. For example, in QCD 
factorization, there are large theoretical uncertainties related to the modelling of power corrections 
corresponding to weak annihilation effects and the chirally-enhanced power corrections to hard 
spectator scattering. Even for two-body B decays, power corrections are of order (10-20)% for 
tree-dominated modes, but they are usually bigger than the central values for penguin-dominated 
decays. Needless to say, 1/nib power corrections for three-body decays may well be larger. However, 
as stressed in Introduction, in this exploratory work we use the phenomenological factorization 
model rather than in the established theories based on a heavy quark expansion. Consequently, 
uncertainties due to power corrections, at this stage, are not included in our calculations, by 
assumption. In view of such shortcomings we must emphasize that the additional errors due to 
such model dependent assumptions may be sizable. 

From Table IIIII we see that the predicted rates for resonant and nonresonant components are 
consistent with experiment within errors. The nonresonant contribution arises dominantly from 
the transition process (88%) via the scalar-density-induced vacuum to KK transition, namely, 
(K + K~\ss\0) , and slightly from the current-induced process (3%). Therefore, it is natural to 
conjecture that nonresonant decays could also play a prominent role in other penguin dominated 
3-body B decays. 

The K + K~ Ks mode is an admixture of CP-even and CP-odd components. By excluding the 
major CP-odd contribution from 4>Ks, the 3-body K + K~ K$ final state is primarily CP-even. The 
K + K~ mass spectra of the B° — ► K + K~ K$ decay from CP-even and CP-odd contributions are 
shown in Fig. [2 For the CP-even spectrum, there are peaks at the threshold and mx+K- = 1-5 
GeV region. The threshold enhancement arises from the /o(980)iT,s and the nonresonant ftf +K 
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FIG. 2: The K + K~ mass spectra for B° — > K + K~ K$ decay from (a) CP-even and (b) CP-odd 
contributions. The insert in (b) is for the eft region. The full K + K~K$ spectrum, which is the sum 
of CP-even and CP-odd parts, measured by BaBar [3] is depicted in (c). 



contributions [see Eq. (|2.18|) ]. 2 For the CP-odd spectrum, the peak on the lower end corresponds 
to the (j)Ks contribution, which is also shown in the insert. The b — > u transition is governed by 
the current-induced process (B K+K ) x (0 K~) [see Eq. ([£4]) ]. From Eq. ([22]) it is clear 
that the b — » u amplitude prefers a small invariant mass of K + and K° and hence a large invariant 
mass of K + and K~. In contrast, the 6 — > c amplitude prefers a small S23. Consequently, their 
interference is largely suppressed. The full K + K~ K$ spectrum, which is the sum of the CP-even 
and the CP-odd parts, has been measured by BaBar [Fig. 0(c)]. It clearly shows the phenomenon 
of threshold enhancement and the scalar resonances Xo(1550) and XcO- 

The decay B -► K S K S K S is a pure penguin-induced mode [cf. Eq. (lA7h ] and it receives 
intermediate pole contributions only from the iso-singlet scalar mesons such as /n(980). Just like 
other KKK modes, this decay is governed by the nonresonant background dominated by the a NR 
term defined in Eq. (|2.18p . Hence, this mode is ideal for determining the unknown parameter cr NR 
which is given in Eq. (|2.26[) . Time-dependent CP violation in neutral 3-body decay modes with 
fixed CP parity was first discussed by Gershon and Hazumi 5l| . 

Results for the decay rates and CP asymmetries in B° — > K + K~ Kgr^, KgKsKgr^ are dis- 
played in Table HVl and Table N] respectively. (For the decay amplitudes of B° — > KsKsKs(l), see 



40 ] for details.) The mixing- induced CP violations are defined by 

2 / Im(e- 2 ^ A C P±A* CP± )dsM<ls 



Skkk,cp± 
Skkk 



12US23 



J \A C p±\ 2 dsi 2 ds 2 3 + J \A C p±\ 2 ds 12 ds 2 3 ' 
2 / Im(e- 2i/3 AA*)ds 12 ds 2 3 

J \A\ 2 dsi 2 dS23 + J |^4| 2 (isi2dS23 

/+ Skkk,cp+ + (1 - /+) Skkk,cp-, 



(2.27) 



2 In our previous work 40j we have argued that the spectrum should have a peak at the large m K + K ~ 
end. This is because we have introduced an additional nonresonant contribution to the w_ parameter 
parametrized as cj^ r — k 2pj |' P2 and employed the B~ — > D°K°K~ data and applied isospin symmetry 

S 12 

to the B —> KK matrix elements to determine the unknown parameter k. Since this nonresonant term 
favors a small m K + Ks region, a peak of the spectrum at large m K + K - is thus expected. However, such a 
bump is not seen experimentally 16]. In this work we will no longer consider this term. 
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TABLE IV: Branching ratios for B° — > K + K~ K$, KgKgKs, K$KsKl decays and the fraction 
of CP-even contribution to B° — ► K + K~Ks, /+. The branching ratio of CP-odd K + K~K$ with 
excluded is shown in parentheses. Results for (K + K~ Kl)cp± are identical to those for 



(K + K Ks)cp^- For theoretical errors, see Table ITTTl Experimental results are taken from 501 ] . 



Final State i3(10- b ) theory i3(10- b ) expt 



K+K~K S 
(K+K-K s )cp- 


Q on+0. 19+2.28+0.07 
y - oy -0.21-l. 81-0.08 

qo+0.10+1.82+0.05 
°- oo -0. 12-1.49-0.06 

1 C7+0. 09+0. 46+0. 02 
J-- ' -0.10-0.32-0.02 

(Ci 1/1+0.06+0.14+0.01 \ 
v u — 0.06— 0.06— 0.01 / 


12.4 ±1.2 


KsK s K s 


input 


6.2 ±0.9 


K S K S K L 


7 «Q+0. 01+1. 37+0.03 
' - uo -0. 01-1. 19-0.03 


< 14 




^theory 
J + 


i-expt 
J + 


K+K'Ks 


H qq+0. 01+0. 01+0. 00 
u - yo -0. 01-0. 02-0. 00 


0.91 ±0.07 


^•theory 


K + K-K L 


n qo+0. oi+o. oi+o. oo 

u - yo -0. 01-0. 02-0. 00 





where A is the decay amplitude of B° — » K + K Kgr^ or KsKsK^^l) an d A is the conjugated B° 
decay amplitude, and /+ is the CP even fraction defined by 



r + r 



(2.28) 

excluded. 



Generally, it is more convenient to define an effective sin 2(3 via Sf = —77/ sin 2/3 e g- with 77/ = 2/+ — 1 
for K + K~ Ks- The predicted value of /+ is consistent with the data but it is on the higher end of 
the experimental measurement because the CP-odd contributions from the vector mesons p, oj, ■ ■ ■ , 
are OZI suppressed and the CP-odd nonresonant contribution is constrained by the 7r + ir~ir~ rate. 
The deviation of the mixing-induced CP asymmetry in B° — > K + K~ K$ and KsK$Ks from 



that measured in B — > cg ifg, i.e. sin2/?0 c _^ s = 0.681 ±0.025 [50J, namely, Asin2/? e g- = sin2/3 c g — 
sin 2/3^,-ifg , is calculated from Table M to be 

ABm20 K+K - Kg = 0.047i°;°|, 

Asm2f3 KsKsKs = 0.038l°;°i- (2.29) 

The corresponding experimental values are 0.049 ± 0.10 and —0.101 ± 0.20, respectively. Due to 
the presence of color-allowed tree contributions in B° — > K + K~ K$, it is naively expected that this 
penguin-dominated mode is subject to a potentially significant tree pollution and hence A sin 2/3 e g- 
can be as large as O(10%). However, our calculation indicates the deviation of the mixing-induced 
CP asymmetry in B° — > K + K~ K$ from that measured in P° — > <ft C cKs is very similar to that of 
the KsKgKs mode as the tree pollution effect in the former is somewhat washed out. Nevertheless, 
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TABLE V: Mixing-induced and direct CP asymmetries sin2/3 e ff (top) and Af (in %, bottom), 
respectively, in B° — > K + K~K$ and KsKsKs decays. Experimental results for K + K~Ks and 
K + K~ Kl modes are obtained from the data of B° — > K + K~ K° . Results for (K + K~ Kl)cp± are 
identical to those for (K + K~ Kg)cp^- F° r theoretical errors, see Table Hill Experimental results 
are taken from 



5C|- 



Final state sin2/3 e ff Expt. 



(K + K Ks)<f>K s excluded 


n 7OS+0-001+0.002+0.009 
U. '^o_ 002-0 .001-0.020 


0.73 ±0.10 


(K + K-K s )cp+ 


n 7QQ+0.003+0.006+0.009 
U. rOZ_ 004-0.004-0.020 




(K+K^K^^Kl excluded 


n 7OO+0.001+0.002+0.009 
U. 1 ^o_ .002-0.001-0.020 


0.73 ±0.10 


KsKsKs 


n 71 Q+0.000+0.000+0.008 
u - ' iy -0.000-0.000-0.019 


0.58 ±0.20 


K S K S K L 


n 71 O+0.000+0.000+0.008 
u - ' io -0.000-0. 000-0.019 





A f (%) Expt. 



(K + K Ks)<j>Ks excluded -4.631 l!ol-o!54-o!34 -7 ± 8 

(K^K K S )CP+ -4.86_ 1 .09-0.55-0.35 

{K + K~Kl)^>K l excluded —4.631 \ oilo 54^0 34 - 7 ± 8 

KsKsKs 0.69inilH^ni 14 ± 15 

KsK s K L 0.77lHilH^:^ 



direct CP asymmetry of the former, being of order —4%, is more prominent than the latter. 3 



B. B~ -» ifJTiT decays 

The B~ — > K + K~ K~ decay amplitude has a similar expression as Eq. (|A4|) except that one 
also needs to add the contributions from the interchange S23 — * S12 and put a factor of 1/2 in the 
decay rate to account for the identical particle effect. 

Branching ratios of resonant and nonresonant contributions to B~ — > K + K~ K~ are shown in 
Table PVTl It is clear that the predicted rates of resonant and nonresonant components are consistent 
with the data except for the broad scalar resonance .Xo(1550). Both BaBar and Belle have seen 
a large fraction from X (1550), (121 ± 19 ± 6)% by BaBar [8] and (63.4 ± 6.9)% by Belle Q, 4 
while our prediction is similar to that in B — ► K + K~ K . It is not clear why there is a huge 



In our previous work , A sin 2/3 e ff is found to be 



Asin2/W-ir s = 0.061^4 , Aam2f3 KsKsKs = 0.06lg;{g , 

for sm2j3j/^K s = 0.687 ± 0.032, while direct CP asymmetry is less than 1% in both modes. Note that 
due to an oversight the experimental error bars were not included in our previous paper for the theoretical 
calculation of Asin2/3 e g. 

Belle Q actually found two solutions for the fraction of X (1550)K-: (63.4 ± 6.9)% and (8.21 ± 1.94)%. 
The first solution is preferred by Belle. 
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disparity between B~ — > K + K~ K~ and -E?° — > i ; T + i ; r~A ;0 as far as the Xo(1550) contribution is 
concerned. Obviously, a refined measurement of the Xo(1550) contribution to the K + K~ K~ mode 
is urgently needed in order to clarify this issue. Our result for the nonresonant contribution is in 
good agreement with Belle, but disagrees with BaBar. Notice that Belle did not see the scalar 
resonance /o(980) as Belle employed the E791 result 53] for gf°^ KK which is smaller than g*°~ >7r7r . 
In contrast to E791, the ratio gf°^ KK / gfo^^ 71 is measured to be larger than 4 in the existing e + e~ 
experiments 0j S] 



TABLE VI: Branching ratios (in units of 10 6 ) of resonant and nonresonant (NR) contributions to 
B~ K + K~K~. For theoretical errors, see Table ITTTl 



Decay mode BaBar ]8] Belle [9] Theory 



cj)K- 


4.14 ±0.32 ±0.33 


4.72 ± 0.45 ± 0.35l[];| 


9 q+0.0+0. 5+0.0 
^• y -0.0-0.5-0.0 


/o(980)K- 


6.5 ± 2.5 ± 1.6 


< 2.9 


7 n+0.0+0.4+0.1 
'• u -0.0-0. 7-0.1 


X (1550)K~ 


43 ± 6 ± 3 




1 -1+0.0+0.2+0.0 
J - J --0.0-0.2-0.0 


fo(1710)K- 


1.7 ± 1.0 ±0.3 






NR 


50 ± 6 ± 4 


24.0 ± 1.5 ± 1.8±£;? 


9 r o+0. 9+4.8+0.3 
zo - -l. 0-4.4-0.3 



Total 35.2 ±0.9 ± 1.6 32.1 ±1.3 ±2.4 25.51 



We next turn to the decay B —* K K$Ks- Following 5a], let us consider the symmetric state 
of K°K° 

\K°K°) sym = [\K\ Pl )K\ P2 )) + \K°( Pl )K (p 2 ))] /V2 

= [\K s (pi)K s (p2)) ~ \K L ( Pl )K L (p 2 )}} /V2. (2.30) 

Hence, 

B(B~ K-KsKs) = \[B{B~ -► K-K S K S ) + B(B~ -» K-R L K L )] 

= ~B(B- -^K-(K K°) sym ). (2.31) 

The factorizable amplitude of — > K~K°K° is given by Eq. (|A8j) . Just as other KKK modes, 
this decay is also expected to be dominated by the nonresonant contribution (see Table IVII|) . 
The calculated total rate is in good agreement with experiment. Just as the pure penguin mode 
KsKsKs, the decay B~ — > K~K$Ks also can be used to constrain the nonresonant parameter 

As pointed out in [55J], isospin symmetry implies the relation 

A{B~ -> K-K°K°) = -A(B° -» K°K + K-). (2.32) 

This leads to 

B(B~ -> K-(K°K°) sym ) = T¥L-Lb(B* K + K-k\ k excluded . (2.33) 

Experimentally, this relation is well satisfied: LHS=(23.0±2.6) x 10~ 6 and RHS=(22.1±2.1) x 10~ 6 . 
Hence, the isospin relation Eq. (|2.32j) is well respected. 
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TABLE VII: Branching ratios (in units of 10 6 ) of resonant and nonresonant (NR) contributions 
to B~ — > K~K$Ks- For theoretical errors, see Table ITO1 



Decay mode 


/o(980)if- X (1550)K- NR 


total 


Theory 
Expt. 


r 9+0.0+0.3+0.1 n q 9 +0.00+0.16+0.00 1 9 ^+0.2+2.1+0.1 
J - z -0. 0-0.5-0.1 u - yz -0.00-0.15-0.00 1 ^-^-0.3-2. 0-0.1 


-19 9+O.O+I.5+O.O 
lz - z -0.0-l. 7-0.0 

11.5 ± 1.3 



III. B Ktttt DECAYS 

In this section we shall consider five B — > Ktttt decays, namely, B~ — > if~7r + 7r~, if°7r~7r°, 
B — > K~tt + tt , if tt + tt~ and if V°7T . They are dominated by b — > s penguin transition and 
consist of three decay processes: (i) the current-induced process, (B — » 7r7r) x (0 — ► if), (ii) the 
transition processes, {B — > tt) x (0 — » 7rif), and (-B — > if) x (0 — > vr7r), and (iii) the annihilation 
process (2? — > 0) x (0 — > Ktttt). 

The factorizable amplitudes for i? _ — ► if~7r + 7r~, if tt~tt°, B — > if _ 7r + 7r , if 7r + 7r~ and 
if°7r 7r are given in Eqs. (|A10IIA14|) . respectively. All five channels have the three-body matrix 
element {tttt\ (qb)v~A\B) which has the similar expression as Eqs. (12. 3p and (12. 4p except that the 
pole B* is replaced by B* and the kaon is replaced by the pion. However, there are additional 
resonant contributions to this three-body matrix element due to the intermediate vector p and 
scalar /o mesons 

(7T + ( P2 )TT-(p 3 )\(ub) v _ A \B-) R = E ,_^_ + * r Y,e*-ijP2-pM\(ub) v _ A \B-) 



pol 



+ E 2 ^ - r (foi\W) v _ A \B-), (3.1) 

where denote generic p-type vector mesons, e.g. p = p(770), p(1450), p(1700), • • - . Applying Eqs. 
(fBT|) and |B6j we are led to 

(TT + (p 2 )TT-( P3 )\(ub) V - A \B-) R (K-( P1 )\(SU) V _ A \0) 



= ? T - _ — ^ («i2 - -13) [(m fl + m„)Af « (, 
tub + m Pi 

Likewise, the 3-body matrix element (if ~ir + \(sb) v _ A \B ) appearing in B — » if ~tt + tt° also receives 
the following resonant contributions 

(K-( Pl )TT + (p 2 )\(- S b) v „ A \B-y = E^^— r-E £ * • (pi-P2)(K*%b) v _ A \B°), 



mj <: -s 12 -im K: T Kt pQl 



(3.3) 
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with K* = iT*(892), 1^(1410), -FT*(1680),---. 

For the two-body matrix elements (tt + K~\ (sd)v-A\0), (tt + tt~\(uu)v~a\0) and (7r + 7r~|ss|0), we 
note that 

(K-( Pl )ir+(p 2 )\(sd) v _ A \0) = (7r + (p 2 )\(sd) v _ A \K + (- Pl )} = {px-p2) tl F^{s n ) 



+ 



(Pi+P2) M -F^{s 12 )+F^{s l2 ) , (3.4) 



Kn, 



812 



where we have taken into account the sign flip arising from interchanging the operators s <-> d. 
Hence, 

(K-( Pl )7r + (p 2 )\(sd) v ^ A \0)(7T-(p 3 )\(db) v _ A \B-) 

(m| - m\){m\ - m\) 



F^{s l2 )F^(s l2 ) 



•S23 - S13 



S\2 



(m B - m 7T )(m K - mj 



S12 



However, the form factor F\ also receives resonant contributions 

1 



A K*tvK „,2 



■ — — — fnK*fK*-\ o — — 

1 "~ - si2 - im K fT K * 1 1 mjf*. - s 12 - im^XK^ 



(3.5) 



fK* i {Pl-P2)n , (3.6) 



with 



eJA^ = (K-( Pl )7T + ( P2 )\K*) = g K *^ K e* • (pi - pa), (3.7) 
where ifp? = -^o(1430), ■ • -. Hence, the resonant contributions to the form factor F^ are 



m K *f K *g + 



JK*9 °* 



mi-, — s — irriK*^K* m K , — s — ira^* Tk 
In principle, the weak vector form factor i^ 17r+7r defined by 

(vr + (p 7r+ )vr~(p 7r -)|n7 At n|0) = (p w+ - ) At F 7r+7r , 



(3. 



(3.9) 



can be related to the time-like pion electromagnetic form factors. However, unlike the kaon case, 
the time-like e.m. form factors of the pions are not well measured enough allowing us to determine 
the resonant and nonresonant parts. Therefore, we shall only consider the resonant part which has 
the expression 



mi, — s — im n V c 



(3.10) 



Following Eq. (|2.18p . the relevant matrix elements of scalar densities read 



(vr + (p 2 )vr-(p 3 )|s S |0)=^-2 



mj - s 23 - im h .T h . 



+ (7T + ( P2 )7r-(p 3 )\ss\0) NR , (3.11) 



and 



f n KX—>K 7T"*" 

(K-( Pl )n+(p 2 )\sd\0) = 2 0< 0i ■ TT- + (K-( Pl )n + (p 2 )\sd\0) NR . (3.12) 

i m KT- - S 12 ~ m K* n lK; 
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Note that for the scalar meson, the decay constants fs and fs are defined in Eq. (|Bip and they 
are related via Eq. (|B2|) . The nonresonant contribution (tt + (p 2 )n~ (p3)\ss\0) N R vanishes under the 
OZI rule, while under SU(3) symmetry 5 

(K-( Pl )n + (p 2 )\sd\0) NR = (K + ( Pl )K-(p 2 )\ss\0) NR = f™(s 12 ), (3.13) 

with the expression of f^ R given in Eq. (|2.18p . 

It is known that in the narrow width approximation, the 3-body decay rate obeys the factoriza- 
tion relation 

T{B -» RP -> P 1 P 2 P) = T(B -» RP)B{R -» P 1 P 2 ) 1 (3.14) 

with R being a resonance. This means that the amplitudes A{B —* RP — > P\P 2 P) and A(B — > i?P) 
should have the same expressions apart from some factors. Hence, using the known results for 
quasi-two-body decay amplitude A(B — ► -RP), one can have a cross check on the three-body decay 
amplitude of B — > RP — * P\P 2 P. For example, from Eq. (|A12p we obtain the factorizable 
amplitude A(B° -> £T *° (1430)vr° ; (1430) -► ^-vr + ) as 

(K-{ Pl )i: + {p 2 )^(p z )\T p \^) Kt)0{im = 

/.if 5 (m| -mi. )), (3.15) 



3 / 

a2<V + 2^ a9 — ° 7 



where 



The expression inside {• • •} is indeed the amplitude of P° — ► -K"q (1430)7t given in Eq. (A6) of 
48]. 

The strong coupling constants such as gP^ n+7T and gfo( 9SO )^ n+ 71 are determined from the 
measured partial widths through the relations 

t-i P c 2 r ^ Pc 2 fo 1 >7\ 

rS = ta| 9 ™' ry= 3 4^<^^' (3 - 17) 

for scalar and vector mesons, respectively, where p c is the cm. momentum. The numerical results 
are 

gP-****- = g Q) g K*^K+-K- = 4>59) 

5 /o(980)-7r+7r- = L33 +0.29 GeV) ^o*-^ + ^ = 3 . 8 4 GeV. (3.18) 
In determining the coupling of fo — > ir + ir~ , we have used the partial width 

r(/ (980) -» 7T+7T-) = (34.2±g;|±|;|) MeV (3.19) 



5 The matrix elements of scalar densities can be generally decomposed into D-, F- and 5(singlet)-type 
components. Assuming that the singlet component is OZI suppressed, SU(3) symmetry leads to, for 
example, the relation (Kn\sq\0) NR = (KK\ss\0) NR . 
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measured by Belle [561 ] . The momentum dependence of the weak form factor F K7T (q 2 ) is 
parametrized as 

1 - g 2 /A x 2 + i^/A, 



*"*V) = ; 2/A 2 - L r, , (3-20) 



where A x ~ 830 MeV is the chiral-symmetry breaking scale 57J and r# is the width of the relevant 
resonance, which is taken to be 200 MeV [381 ] . 

The results of the calculation are summarized in Tables IVIIItiXIIl We see that except for 
fo(980)K, the predicted rates for K*tt, iTo(1430)7r and pK are smaller than the data. Indeed, the 
predictions based on QCD factorization for these decays are also generally smaller than experiment 
by a factor of 2~5. This will be discussed in more details in Sec. VI. 



TABLE VIII: Branching ratios (in units of 10 6 ) of resonant and nonresonant (NR) contributions 
to B~ — ► K~tt + tt~ . For theoretical errors, see Table HTT1 



Decay mode 


BaBar [6] 


Belle [7] 


Theory 


K TT 

K* °(U30)tt- 
p°K- 
/o(980)if- 
NR 


9.04 ± 0.77 ± 0.53±g;|^ 
34.4 ± 1.7 ±1.8t?:4 
5.08 ± 0.78 ± 0.39±g;|| 
9.30 ± 0.98 ± 0.5lig;^ 
2.87 ± 0.65 ± 0.43ig;i 


6.45 ± 0.43 ± 0.48lg j| 
32.0 ± 1.0 ± 2.4±i:J 
3.89 ± 0.47 ± 0.29lg i 
8.78 ± 0.82 ± 0.65iS;e4 
16.9 ± 1.3 ± 1.3tJ;J 


n+o.o+o.8+o.o 

°- u -0.0-0.7-0.0 
in e+0.0+3. 2+0.0 
iu -°-0.0-2. 7-0.1 

1 o + O.O+l. 9 + 0.1 
i,a -0.0-0.7-0.1 
7 7+0.0+0.4+0.1 
' • ' -0.0-0.8-0.1 

1 a 7+0.5+11.0+0.2 
io -'-0.6- 6.3-0.2 



Total 64.4 ±2.5 ±4.6 48.8 ±1.1 ±3.6 45.0; 



TABLE IX: Same as Table EHD except for the decay B~ -> K tt~tt° . 



Decay mode 


Theory 


Decay mode 


Theory 


K*-tt° 
K*-{U30)tt° 
p K 


-i c+0. 0+0.3+0. 2 

^-o.o-o.a-o^ 

c c+O.0+1.6+0.1 
d - d -0.0-l. 4-0.1 
1 q+0.0+3.0+0.0 

^^-o.o-o.g-o.o 


K TT 

K* °(U30)tt- 
NR 


1 r+0.0+0.4+0.0 
i>o -0.0-0.3-0.0 
r 9+0.0+1.6+0.0 
d -^-0.0-l. 4-0.0 
ln n+0-2+7.1+0.0 
U-U — 0.2— 3.7— 0.0 


Total 


97 n +0. 3+15. 4+0.2 
z '- u -0.2- 8.8-0.2 





While Belle has found a sizable fraction of order (35 ~ 40)% for the nonresonant signal in 
K~TT + TT~ and K°tt + tt~ modes (see Table H|), BaBar reported a small fraction of order 4.5% in 
K~TT + TT~. The huge disparity between BaBar and Belle is ascribed to the different parameteri- 
zations adopted by both groups. BaBar used the LASS parametrization to describe the Ktt 
S-wave and the nonresonant component by a single amplitude suggested by the LASS collaboration 
to describe the scalar amplitude in elastic Kit scattering. As commented in [7J, while this approach 
is experimentally motivated, the use of the LASS parametrization is limited to the elastic region 
of M(Ktt) < 2.0 GeV, and an additional amplitude is still required for a satisfactory description 
of the data. In our calculations we have taken into account the nonresonant contributions to the 
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TABLE X: Same as Table [yTTTl except for the decay B — » K 7r + 7r . 



Decay mode Belle [131 Theory 





5.6 ±0.7 ±0.5lg;| 


9 1+0.0+0.5+0.3 
z - i -0.0-0.5-0.3 


i<r *-(1430> + 


30.8 ± 2.4 ± 2.4±§:g 


in 1+0.0+2.9+0.1 
iu - i -0.0-2.5-0.2 


p°K° 


6.1 ± 1.0±0.5lJ;S 


9 n+0.0+1. 9+0.1 
z - u -0.0-0.9-0.1 


/o(980)K° 


7.6 ± 1.7 ± 0.7±g;f 


7 7+0.0+0.4+0.0 
'•'-0.0-0.7-0.0 


NR 


19.9 ± 2.5 ± 1.6l?;^ 


i rr fi +0.1+8. 3+0.0 
io -°-0.1-4.9-0.0 



Total 47.5 ± 2.4 ± 3.7 42.01; 



TABLE XI: Branching ratios (in units of 10 ) of resonant and nonresonant (NR) contributions 



to B° — ► K 7T + 7r°. Note that the branching ratios for K* tt + and i^*°7r° given in [14] and 



151 ] are their absolute ones. We have converted them into the product branching ratios, namely, 



B(B —* Rh) x B(R — > hh). For theoretical errors, see Table HTT1 



Decay mode 


BaBar [14] 


Belle [15J 


Theory 


K*~TT + 


3.6 ±0.8 ±0.5 


a q+1. 5+0. 5+0.8 
^• y -l. 5-0. 3-0.3 


1 n+o.o+o.3+o.i 
i - u -o. 0-0.3-0.1 


ifV 


2.0 ±0.6 ±0.3 


< 2.3 


1 n+0-0+0.3+0.2 

J-u-o.o-o.a-o.i 


i^o"(1430)^+ 


11.2 ±1.5 ±3.5 


5.i ± i.5ig:? 


c n+0.0+1. 5+0.1 
°- u -o. 0-1. 3-0.1 


K*°(1430)7r° 


7.9 ±1.5 ±2.7 


c 1+1.6+0.5 


a 9+O.O+I.4+O.O 
*- z -0. 0-1. 2-0.0 


p+K- 
NR 


8.6 ±1.4 ±1.0 
< 4.6 


i e 1+3.4+1.4+2.0 
iJ - L -3.3-1.5-2.1 
e 7 +2.7+0.5 .04 
°- ' -2.5-0.4 ^ y - 4 


9 r+0.0+3.6+0.2 
z - d -0.0-1.4-0.2 
n e+0.3+6.6+0.0 
y - D -0. 2-3.5-0.0 


Total 


34.9 ±2.1 ±3.9 


36.6^1 ± 3.0 


90 n+0.2+16.1+0.2 
zo.y_ Q 2 _ g 4 _o 2 



two-body matrix elements of scalar densities, {K7r\sq\0) . Recall that a large nonresonant contribu- 
tion from (KK\ss\0) is needed in order to explain the observed decay rates of B° — > KgKgKg and 
5~ — > K~ K$K$- From Tables IVIimXIII we see that our predicted nonresonant rates are in agree- 
ment with the Belle measurements. The reason why the nonresonant fraction is as large as 90% in 
KKK decays, but becomes only (35 ~ 40)% in Kirn channels (see Table IJ) can be explained as 
follows. Under SU(3) flavor symmetry, we have the relation (Kir\sq\0) NR = (KK\ss\0) NR . Hence, 
the nonresonant rates in the K~ir + ir~ and K°tt + tt~ modes should be similar to that in K + K~K° 
or K + K~K~ . Since the KKK channel receives resonant contributions only from <f> and fo i mesons, 
while K* , KfojPi, fa resonances contribute to Kirir modes, this explains why the nonresonant frac- 
tion is of order 90% in the former and becomes of order 40% in the latter. Note that the predicted 
nonresonant contribution in the K~ir + ir° mode is larger than the BaBar's upper bound and barely 
consistent with the Belle limit. It is conceivable that the SU(3) breaking effect in (Kir\sq\0) NR 
may lead to a result consistent with the Belle limit. 

It is interesting to notice that, based on a simple fragmentation model and SU(3) symmetry, 
Gronau and Rosner 



551 ] found the relations 



T(B- -> K + K-K-) NR = 2T(B° -> K + K~ K°) NR = 2T(B~ -» JTtt+tOnr 
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TABLE XII: Same as Table Em] except for the decay B -> K\ tt . 



Decay mode 


/ (980)K U A*V Kq U (1430)^° NR 


Total 


Theory 


q c+0.0+2. 0+0.0 n cc;+0.00+0.16+0.00 9 q+0.0+0. 8+0.0 c q+0.0+1. 8+0.0 
°-°-0.0-0. 4-0.0 u - d,J -0.00-0.13-0.00 z -°-0.0-0. 6-0.0 d -°-0.0-l. 1-0.0 


1 9 Q+0.0+4. 0+0.1 

lz - y -0.0-3.0-0.1 



TABLE XIII: Branching ratios, mixing-induced and direct CP asymmetries for B — ► i^57r + 7r 
decays. Results for {Kl-kit)cp± are identical to those for (AT,s7r7r)cp=p For theoretical errors, see 
Table HH 



Final state 


Branching ratio 


(K S K + TT~)CP+ 

(K s ir + ir~)cp- 
/+ 


1 3 =9+0.02+4.03+0.01 
io - JZ -0.03-3.06-0.01 
7 ^+0.10+3.79+0.02 
'•^°-0. 08-2. 32-0.02 
n ctx+0. 00+0. 03+0. 00 
U - DO -0. 00-0. 04-0. 00 


Final state 


sin 2/?eff 


(K s -K + ir~)cp+ 
(K57r + vr-) ful i 


n RQQ+O.000+0. 003+0.003 

u.oyo_ .000-0 .002-0 .014 

n 71 O+0.001+0. 017+0.008 
u - ' lo -0.001-0. 007-0.018 


Final state 


A f (%) 


(K s n + 7r~)cp+ 
(K 5 vr+7r-) ful i 


a 97+0.00+0.19+0.28 
^• z ' -0.00-0.12-0.35 
A QA +0.03+0.03+0.32 
*- y ^-0. 02-0. 05-0. 40 



2r(£?° -> aV+Onr = 4T(B° -» ^-vr + 7r°) NR . (3.21) 



Again, a large nonresonant background in A~7r + 7r _ and A'°7r + 7r~ is favored by this model. 

Although the B° — > As7r 7r rate has not been measured, its time-dependent CP asymmetries 
have been studied by BaBar [58|] with the results 

sin 2/3 eS = -0.72 ±0.71 ±0.08, A C p = -0.23 ± 0.52 ± 0.13 . (3.22) 

Note that this mode is a CP-even eigenstate. We found that its branching ratio is not so small, 
of order 6 x 10 -6 , in spite of the presence of two neutral pions in the final state (see Table IXlT)) . 
Theoretically, we obtain 

cl-nOR — n 79Q+0.000+0.001+0.009 a _ f n 90+0. 09+0. 07+0.02 \ 0/ fo oq\ 

sm2p eS - 0.729„ 000-0.001-0.020' a cp - ^0.28_ .o6-o.o6-o.o2 ) /o- K 6 - 26 ) 

Finally, we consider the mode Kstt + tt" which is an admixture of CP-even and CP-odd compo- 
nents. Results for the decay rates and CP asymmetries are displayed in Table IXIIII We see that 
the effective sin 2/3 is of order 0.718 and direct CP asymmetry of order 4.9% for Ks^ + tt^ ■ 



IV. B -> KKir DECAYS 



We now turn to the three-body decay modes dominated by b — ► u tree and b — * d penguin 
transitions, namely, KKtt and tttttt. We first consider the decay B~ — > K + K~ir~ whose factorizable 
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TABLE XIV: Same as Table IVlTIl except for the decay B 



K+K--K- 



Decay mode 


/o(980)vr- K*°K~ Kjp(1430)K- NR 


Total 


Theory 
Expt. 


n cn+0-00+0.06+0.02 n 9 o+0.00+0.04+0.02 n no+0.00+0. 18+0.09 i H +0. 5+0.4+0. 2 
U - OU -0.00-0.04-0.02 u - ZJ -0.00-0.04-0.02 U - OZ -0.00-0.16-0.08 L -°-0. 5-0.2-0. 2 


A n+0.5+0.7+0.3 
^• u -0. 6-0.5-0. 3 

< 6.3 (BaBar) [59] 

< 13 (Belle) [U] 



TABLE XV: Same as Table IVTnl except for B — > ir + 7r ir . The nonresonant background is used 
as an input to fix the parameter a NR defined in Eq. f|2.8[) . 



Decay mode 


BaBar [5] 


Theory 


p°7T- 


8.8 ±1.0 ±0.6±°/ 7 


7 7+0.0+1.7+0.3 
' -'-0.0-1.6-0.2 


/o(980)vr- 


1.2 ± 0.6 ± 0.1 ± 0.4 < 3.0 


n qq+o. oo+o. 01+0.03 

u - oy -0. 00-0. 01-0. 02 


NR 


2.3 ±0.9 ±0.3 ±0.4 < 4.6 


input 


Total 


16.2 ± 1.2 ±0.9 


19 n +l. 1+2. 0+0.4 
iz - u -l. 2-1. 8-0.3 



amplitude is given by Eq. (1A9I) . Note that we have included the matrix element {K + K~\dd\0). 
Although its nonresonant contribution vanishes as K + and K~ do not contain the valence d or d 
quark, this matrix element does receive a contribution from the scalar /o pole 

(K+(p 2 )K-( P3 )\dd\0) R = J2 2 /0 ' /0 ' ■ 7— . (4.1) 

i m f 0l ~ s 23 - irn fo T fo . 

where (/ o |dd|0) = rnf Q ff Q . In the 2-quark model for / (980), // ( 980 ) = // (980) sin6>/^2. Also note 
that the matrix element (K~ {p 3 )\{sb)v -a\B^) (ir~ {p\)K + (p 2 )\(ds)v~A\0) has a similar expression 
as Eq. (13. 5j) except for a sign difference 

(K-(p 3 )\(sb)v~A\B-)(TT-( Pl )K+(p 2 )\(ds)v-A\0) 

(m 2 B - m 2 K )(m 2 K - ml) 



F 1 BK (s 12 )F 1 K ^s 12 ) 



S23 ~ Sl3 



S12 

-F BK (s 12 )F«^s 12 ) ^-< ){ <- m *\ (4.2) 

Sl2 



As in Eq. f|3.8[) . the form factor F^ n receives a resonant contribution for the K* pole. 

The nonresonant and various resonant contributions to B~~ — > /C + i^ _ 7r _ are shown in Table 
IXVII The predicted total rate is consistent with upper limits set by BaBar and Belle. 

V. B -> Trvrvr DECAYS 



The factorizable amplitudes of the tree-dominated decay B — > 7r + 7r it and S° — > it + it it 
are given by Eqs. (|A15j) and (|A16h . respectively. We see that the former is dominated by the p° 
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TABLE XVI: Same as Table Era] except for the decay B° -» tt + tt~tt°. 



Decay mode 




p TT + 


p°7T° 


/0(980)7T° 


NR 


Total 


Theory 


a c+0.0+1. 1+0.2 
°- -o. 0-1. 0-0.1 


1 k cr+0.0+4.0+0.3 
±J - J -0.0-3.5-0.3 


i n +0. 0+0. 3+0.0 
1 -0. 0-0. 2-O.0 


n nin+ - 000 +°- 003 + - 000 

U.U1U_ QOO-0 .002-0 .000 


n n^-^+o.oi+o.oo 

u - uo _0. 02— 0.01-0.00 


9fi O+0.0+5.6+0.2 

ZD - -0.0-5.0-0.2 



pole, while the latter receives and p° contributions. As a consequence, the 7r + 7r _ 7r° mode has a 
rate larger than 7r + 7r _ 7r~ even though the former involves a 7r° in the final state. 

The 7r + 7r~7r~ mode receives nonresonant contributions mostly from the b — > it transition as 
the nonresonant contribution in the matrix element (7r + 7r~|d(f|0) is suppressed by the smallness of 
penguin Wilson coefficients a% and as- Therefore, the measurement of the nonresonant contribution 
in this decay can be used to constrain the nonresonant parameter a NR in Eq. (|2.8p 



VI. DIRECT CP ASYMMETRIES 



Direct CP asymmetries for various charmless three-body B decays are collected in Table IXVII1 
Mixing-induced and direct CP asymmetries in B° — ► K + K~ K$,l and KsKsKs : l decays are al- 
ready shown in Table IVl It appears that direct CP violation is sizable in K + K~ K~ and K + K~ir~ 
modes. 

The major uncertainty with direct CP violation comes from the strong phases which are needed 
to induce partial rate CP asymmetries. In this work, the strong phases arise from the effective 
Wilson coefficients a? listed in (|A3p and from the Breit-Wigner formalism for resonances. Since 
direct CP violation in charmless two-body B decays can be significantly affected by final-state 
rescattering [60(|, it is natural to extend the study of final-state rescattering effects to the case of 
three-body B decays. We will leave this to a future investigation. 



TABLE XVII: Direct CP asymmetries (in %) for various charmless three-body B decays. For 
theoretical errors, see Table Hill Experimental results are taken from 5C|. 



Final state 


BaBar 


Belle 


Theory 


K+K-R- 


-2 ±3 ±2 




in A +1-7+0.9+0.9 
lu -^-l. 3-1.0-0.8 


K-K S K S 


-4 ±11 ±2 




q Q + 0. + 0.6 + 0. 3 
a - y -0.0-0.8-0.3 


K+K-tt- 


0±10±3 




17 c+1.9+2.2+0.0 
1 ' - d -3.8-3.4-0.2 


K~TT + TT~ 


-1.3 ±3.7 ±1.1 


4.9 ±2.6 ±2.0 


O O+O. 7 + 0.4+0. 3 

°-°-Q.5-0.4-0.2 


K-TT+TT 




7± 11 ± 1 


c q+0.6+1.4+0.5 
u -°-0. 7-1.4-0. 5 


17° + - 






A q+0. 0+0.0+0.3 
^• y -0. 0-0.1-0. 4 




-23 ± 52 ± 13 


-17 ±24 ±6 


n O o+0.09+0.07+0.02 
U - ZO -0.06-0.06-0.02 


K°ir-ir° 






n a +0.0+0.4+0.0 

u '^-0. 0-0.4-0.0 




-1±8±3 




a a+0. 8+1. 2+0.0 
^-0.6-0.9-0.2 








q n +0. 1+0.2+0.3 
°- u -0. 1-0.3-0.2 
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VII. TWO-BODY B -^VP AND B —> SP DECAYS 

Thus far we have considered the branching ratio products B(B — ► Rh\)B(R — ► ^2^3) with 
the resonance R being a vector meson or a scalar meson. Using the experimental information on 
B(R - ^3) B 

£(A*° -» A+vr") = B{K* + -» A tt + ) = 2£(A*+ -» A+vr ) = -, 

3 

B(iC^ (1430) -> A+tT) = 2£(A*+(1430) -> A+vr ) = -(0.93 ± 0.10), 

3 

B{4>^K + K-) = 0.492 ±0.006. (7.1) 

one can extract the branching ratios of S — » VP and S — ► SP. The results are summarized in 
Table IXVrnl 

Two remarks about the experimental branching ratios are in order: (i) The BaBar results for the 
branching ratios of W -» K*-ir+, K* V, A*-(1430)vr+ are inferred from the three-body decays 
B° — > A°7r + 7r~ (see Table 1X1)) and Belle results are taken from B — ► A~7r + 7r° (see Table |X|) . (ii) 
Branching ratios of B° — > <^A shown in Table IXVIIII are not inferred from the Dalitz plot analysis 
of B — > AAA decays. 

For comparison, the predictions of the QCD factorization approach for B — > UP [61] and 
— > [48j are also exhibited in Table IXVIIII In order to compare theory with experiment for 
B — > /o(980)A — ► 7r + 7r~ A, we need an input for 0(/o(98O) — > 7r + 7r~). To do this, we shall use the 
BES measurement 



45] 



r(A(M0)-„) . 7s «.n. (7 . 2) 



r(/o(980) -. rr) + r(/o(980) - A' A') 



Assuming that the dominance of the /o(980) width by tttt and A A and applying isospin relation, 
we obtain 

B(/ o (980) -» tt+tT) = 0.50lg;S , B(/o(980) -» K + K~) = 0.125ig;{£! ■ (7.3) 
At first sight, it appears that the ratio defined by 

S(/o(980) -. jr+ir") 1 ' 

is not consistent with the value of 0.69 ± 0.32 inferred from the BaBar data (see Tables IVII and 

R _ T{B~ ^/o(980)A-;/ (980) -> K+K~) _ 6.5 ± 2.5 ± 1.6 

r(B~ -» / (980)A-; / (980) -» vr+vr-) 9.3 ± 1.0±°;| ' 1 ' j 

where we have applied the narrow width approximation Eq. (I3.14p . 

The above-mentioned discrepancy can be resolved by noting that the factorization relation Eq. 
(|3.14p for the resonant three-body decay is applicable only when the two-body decays B — ► RP 
and R — > P\Pi are kinematically allowed and the resonance is narrow, the so-called narrow width 
approximation. However, as the decay /o(980) — * K + K~ is kinematically barely or even not 
allowed, the off resonance peak effect of the intermediate resonant state will become important. 
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TABLE XVIII: Branching ratios of quasi-two-body decays B — > VP and B — > SP obtained from the 
studies of three-body decays based on the factorization approach. Unless specified, the experimental 
results are obtained from the 3-body Dalitz plot analyses given in previous Tables. Theoretical 
uncertainties have been added in quadrature. QCD factorization (QCDF) predictions taken from 
61] for VP modes and from [48] for SP channels are shown here for comparison. 



Decay mode 




BaBar 


Belle 


QCDF 


This work 


4>K° 




8 - 4 — i:i ±0 - 5 a 


9.0l?j±0.7 b 


A 1+0.4+1.7+1.8+10.6 
*' J — 0.4-1.6-1.9- 3.0 


K o+l.O 
°-°-0.9 


4>K~ 




8.4 ±0.7 ±0.7 


9.60 ± 0.92±J;g| 


a r+0.5+1. 8+1.9+11.8 
*'°-0.4-1.7-2.1- 3.3 


s.9±i:S 


K*\- 




13.5 ± 1.2±g;| 


9.8 ±0.9±\l 


q c+0.4+1. 5+1. 2+7.7 
J - D -0.3-l. 4-1.2-2.3 


4.4ll; 


K*\° 




3.0 ±0.9 ±0.5 


< 3.5 


n 7 +0.1+0.5+0.3+2.6 
u - ' -0.1-0.4-0.3-0.5 


i r+0.5 
i -°-0.4 


K*--k+ 




11.0 ±1.5 ±0.7 


8-4 ±1.118:1 


q q + 1.4+1. 3 + 0.8 + 6. 2 

o -°-1.2-1.2-0.8-1.6 


q 1+0.9 
°- J — 0.9 


K*-TT° 




6.9 ± 2.0 ± 1.3 b 




O q + 1. 1 + 1. + 0.6 + 4. 4 

o -°-1.0-0.9-0.6-1.4 


9 9+O.6 
z ' z -0.5 


K*°K~ 








n q n +0.11+0.12+0.09+0.57 
U - OU -0.09-0.10-0.09-0.19 


S5+ - 06 
u ' oo -0.06 


p°K- 




5.1 ±0.81™ 


3.89 ± 0A7t° H 


9 R +0.9+3. 1+0.8+4. 3 
z - D -0.9-l. 4-0.6-1. 2 


i o+l.O 
L - 6 -0.7 


p°K° 




4.9 ±0.8 ±0.9 


6.1 ±1.0 ±1.1 


a e+0.5+4. 0+0.7+6.1 
^• D -0.5-2. 1-0.7-2.1 


9 n+1.9 
z - u -0.9 


p + K- 




8.6 ±1.4 ±1.0 


i c 1+3.4+2.4 
10 "- L -3.3-2.6 


7 ,.+1.8+7.1+1.2+10.7 
'•^-1.9-3.6-1.1- 3.5 


9 rr+3.6 
Z -°-1.4 


p K 




8.0±J;| ± 0.5 b 




c O+0.6+7.0+1. 5+10.3 
°-°-0.6-3.3-1.3- 3.2 


1 q + 3.0 


p°TT- 




8.8 db l-0±g;3 


8 -°— i:B ± °- 7 b 


i i q+6.3+3.6+2.5+1.3 
ii - y -5.0-3. l-l. 2-1.1 


7.7+_\l 


p~TT + 








91 9 +10.3+8.7+1.3+2.0 
8.4-7.2-2.3-1.6 


15.5118 


p + TT~ 








i K 4 +8.0+5.5+0.7+1.9 
i0 -^-6.4-4.7-1.3-1.3 


8 - 5 -L0 


A° 




1.4 ±0.6 ±0.3 


o 1+0.9+0.6 

•J^-o.s-o.s 


n A +0.2+0.2+0.9+0.5 
u -^-0.2-0. 1-0.3-0.3 


i n +0.3 
••- u -0.2 


/o(980)K°;/ - 


-* 7T + 7r~ 


5.5 ±0.7 ±0.6 


7.6 ± 1.7lg;g 


a 7+0.1+2.1+2.3 c 
°-' -0.1-1.5-1.1 


7 7+0.4 
'•'-0.7 


/o(980)if-;/ - 


-* ir + ir~ 


9.3 ± 1.0±gJ 


8.8 ±0.8l5:| 


7 O+0.2+2.3+2.7 c 
' -°-0. 2-1.6-1. 2 


7 7+0.4 
'•'-0.8 


/o(980)K°;/ - 


->• K + K~ 


5.3 ± 2.2 






5-8l8:s 


/o(980)if-;/ - 


-» K+K- 


6.5 ±2.5 ±1.6 


< 2.9 




7 .oi8: 7 


/o(980)vr-;/ - 


-> ir + ir~ 


< 3.0 




n r+o. o+o. 2+0.1 c 

u -°-0. 0-0.1-0.0 


o.39i8:8^ 


/o(980)vr-;/ - 


-> K + K~ 








o ^n+ 06 

u - ou -0.04 


/o(980)7r°;/ - 


► 7T + 7T~ 






n n9 +0.01+0.02+0.04 c 
u - uz -0.01-0.00-0.01 


nin+ 0003 

U.U1U_Q QQ2 


Zq (1430)^- 




36.6 ± 1.8 ±4.7 


51.6 ± 1.7l™ 


i i n +10. 3+7.5+49. 9 
ii - u - 6.0-3.5-10.1 


i6.9±|:l 


Zo°(1430)^° 




12.7 ±2.4 ±4.4 


9.8 ±2.5 ±0.9 


fi ,+5.4+2.2+26.1 
°- 4 -3.3-2.1- 5.7 


6-8±?:i 


iT^(1430)7r + 




36.1 ±4.8 ±11.3 


49.7 ± 3.8±|;? 


1 1 q+9. 4+3. 7+45.8 
iJ --°-5.8-3.7- 9.9 


16.2US 


K^-(U30)tt° 








e q+4. 7+1.6+22. 3 
°-°-2.8-1.7- 4.7 


8.9±|:i 


K$°(U30)K- 




< 2.2 6 






i q+0.3 
i - J -0.3 



"From the Dalitz plot analysis of B° — > K + K~K° decay measured by BaBar (see Table ITTT1) . we obtain 

B(B° — > °) = (6.2 ± 0.9) x 1CP 6 . The experimental value of BaBar cited in the Table is obtained from a 

direct measurement of B° — > 
b not determined directly from the Dalitz plot analysis of three-body decays. 
c We have assumed B(/ (980) -> tt+tt") = 0.50 for the QCDF calculation. 
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Therefore, it is necessary to take into account the finite width effect of the /o(980) which has a 
width of order 40-100 MeV Q]. In short, one cannot determine the ratio R by applying the narrow 
width approximation to the three-body decays. That is, one should employ the decays B — * Kirir 
rather than B — > KKK to extract the experimental branching ratio for B — > fo(980)K provided 
B(/o(980) — > 7T7r) is available. 

We now c omp are the present work for B — ► VP and B — ► SP with the approach of QCD 
factorization 34], |48J. In this work, our calculation of 3-body B decays is similar to the simple 



generalized factorization approach 62], |63J by assuming a set of universal and process independent 
effective Wilson coefficients with p = u,c in Eq. (|A3|) . In QCDF, the calculation of a? is 
rather sophisticated. They are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex corrections and hard spectator interactions. In general, 
they have the expressions 34j, l6l| 



a p i (M 1 M 2 



Ci + 



Cj±l 



Ni(M 2 ) + 



Cj±l Cfoi s 



Vi(M 2 ) + —Hi{MxM 2 )\ + P[(M 2 



(7.6) 



where % = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), Cj are the Wilson coefficients, 
C F = (Nl - 1)/(2N C ) with N c = 3, M 2 is the emitted meson and Mi shares the same spectator 
quark with the B meson. The quantities Vi(M 2 ) account for vertex corrections, Hi{M\M 2 ) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 
quark of the B meson and Pi(M 2 ) for penguin contractions. Hence, the effective Wilson coefficients 
a P (M\M 2 ) depend on the nature of Mi and M 2 ; that is, they are process dependent. Moreover, they 
depend on the order of the argument, namely, a p {M 2 M\) ^ a P {M\M 2 ) in general. In the above 
equation, Ni{M 2 ) vanishes for i = 6,8 and M 2 = V, and equals to unity otherwise. For three- 
body decays, in principle one should also compute the vertex, gluon and hard spectator-interaction 
corrections. Of course, these corrections for the three-body case will be more complicated than the 
two-body decay one. One possible improvement of the present work is to utilize the QCDF results 
for the effective parameters a?(MiM 2 ) in the vicinity of the resonance region. 

We next proceed to the comparison of numerical results. For cpK, K*ir and K*K modes, the 
QCDF and the present work have similar predictions. For the p meson in the final states, QCDF 
predicts slightly small pK and too large pir compared to experiment. 6 In contrast, in the present 
work we obtain reasonable pir but too small pK. This is ascribed to the form factor Aq P (0) = 



0.37 ± 0.06 employed in 61] that is too large compared to ours A (0) = 0.28 ± 0.03 (see Table 



IXIXp . Recall that the recent QCD sum rule calculation also yields a smaller one Aq P (0) = 0.30l ' ' 



0.03 



For B -> f (980)K and B -> / (980)vr, QCDF [48(] and this work are in agreement with 
experiment. The large rate of the /o(980)i^ mode is ascribed to the large /o(980) decay constant, 
//o(980) ~ 460 MeV at the renormalization scale p = 2.1 GeV [48] . In contrast, the predicted 
K (1430)tt~ and Kq~(U30)tt + are too small compared to the data. The fact that QCDF leads to 
too small rates for cf)K, K*tt, pK and i^Q(1430)7r may imply the importance of power corrections 



6 Recall that the world average of the branching ratio of B° 
predicts it to be ~ 36.6 x 1CT 6 



61|. 



p ± ir^ is (24.0±2.5) x 1CT 6 [50(, while QCDF 
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due to the non-vanishing pa and pu parameters arising from weak annihilation and hard spectator 
interactions, respectively, which are used to parametrize the endpoint divergences, or due to possible 
final-state rescattering effects from charm intermediate states 60]. However, this is beyond the 
scope of the present work. 



VIII. CONCLUSIONS 

In this work, an exploratory study of charmless 3-body decays of B mesons is presented using 
a simple model based on the framework of the factorization approach. The 3-body decay process 
consists of resonant contributions and the nonresonant signal. Since factorization has not been 
proved for three-body B decays, we shall work in the phenomenological factorization model rather 
than in the established theories such as QCD factorization. That is, we start with the simple idea 
of factorization and see if it works for three-body decays. Our main results are as follows: 

• If heavy meson chiral perturbation theory (HMChPT) is applied to the three-body matrix 
elements for B — > P1P2 transitions and assumed to be valid over the whole kinematic region, 
then the predicted decay rates for nonresonant 3-body B decays will be too large and even 
exceed the measured total rate. This can be understood because chiral symmetry has been 
applied beyond its region of validity. We assume the momentum dependence of nonresonant 
amplitudes in the exponential form e~ a ^ PB '^ Pi+p ^ so that the HMChPT results are recovered 
in the soft meson limit pi, pj — ► 0. The parameter a NR can be fixed from the tree-dominated 
decay B~ — > Tr + ir~ir~. 

• Besides the nonresonant contributions arising from B —* P1-P2 transitions, we have identified 
another large source of the nonresonant background in the matrix elements of scalar densi- 
ties, e.g. (KK\ss\0) which can be constrained from the K$K$Ks (or K~ K$Ks) mode in 
conjunction with the mass spectrum in the decay B~° — > K + K~lC° . 

• All KKK modes are dominated by the nonresonant background. The predicted branching 
ratios of K + K^K S ^, K + K~ K~ and K~K$Ks modes are consistent with the data within 
the theoretical and experimental errors. 

• Although the penguin-dominated B° — ► K + K~ Kg decay is subject to a potentially signifi- 
cant tree pollution, its effective sin 2(3 is very similar to that of the KsK$Ks mode. However, 
direct CP asymmetry of the former, being of order —4%, is more prominent than the latter, 

• The role played by the unknown scalar resonance Ao(1550) in the decay B~ — > K + K~ K~ 
should be clarified in order to see if it behaves in the same way as in the K + K~7C° mode. 

• Applying SU(3) symmetry to relate the nonresonant component in the matrix element 
(Kir\sq\0) to that in (KK\ss\0), we found sizable nonresonant contributions in K~ir + ir~ 
and K°iT + 7r~ modes, in agreement with the Belle measurements but larger than the BaBar 
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results. In particular, the predicted nonresonant contribution in the i^~7r + 7r mode is consis- 
tent with the Belle limit and larger than the BaBar's upper bound. It will be interesting to 
have a refined measurement of the nonresonant contribution to this mode to test our model. 

The 7r + 7T~7r° mode is predicted to have a rate larger than 7r + 7r~7r~ even though the former 
involves a 7r° in the final state. This is because the latter is dominated by the p° pole, while 
the former receives /) ± and p° resonant contributions. 

Among the 3-body decays we have studied, the decay B~ — > K + K~n~ dominated by b —>■ u 
tree transition and b — > d penguin transition has the smallest branching ratio of order 4 x 1CT 6 . 
It is consistent with the current bound set by BaBar and Belle. 

Decay rates and time-dependent CP asymmetries in the decays i^5-7r°7r°, a purely CP-even 
state, and Ks^ + ^, an admixture of CP-even and CP-odd components, are studied. The 
corresponding mixing-induced CP violation is found to be of order 0.729 and 0.718, respec- 
tively. 

Since the decay /o(980) — > K + K~ is kinematically barely or even not allowed, it is crucial 
to take into account the finite width effect of the /o(980) when computing the decay B — > 
fo(980)K — ► KKK. Consequently, one should employ the Dalitz plot analysis of Ktttt mode 
to extract the experimental branching ratio for B — > fo(980)K provided H(/o(980) — > tttt) is 
available. The large rate of B — > fo(980)K is ascribed to the large /o(980) decay constant, 
4(980) ~ 460 MeV. 

The intermediate vector meson contributions to 3-body decays e.g. p, cf>, K* 
are identified through the vector current, while the scalar meson resonances e.g. 
/o(980), Xo(1550), i^o(1430) are mainly associated with the scalar density. Their effects 
are described in terms of the Breit-Wigner formalism. 

Based on the factorization approach, we have computed the resonant contributions to 3- 
body decays and determined the rates for the quasi-two-body decays B — > VP and B — > SP. 
The predicted pir, /o(980)if and /o(980)-7r rates are consistent with experiment, while the 
calculated 4>K, K*ir, pK and -KTq(1430)7t are too small compared to the data. 

Direct CP asymmetries have been computed for the charmless 3-body B decays. We found 
sizable direct CP violation in K + K~K~ and K + K~tt~ modes. 

In this exploratory work we use the phenomenological factorization model rather than in 
the established theories based on a heavy quark expansion. Consequently, we don't have 
l/mft power corrections within this model. However, systematic errors due to such model 
dependent assumptions may be sizable and are not included in the error estimates that we 
give. 
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Note added: After the paper was submitted for publication, BaBar (arXiv:0708.0367 [hep-ex]) has 
reported the observation of the decay B + — > K + K~ir + with the branching ratio (5.0 ± 0.5 ± 0.5) x 
10~ 6 . Our prediction for this mode (see Table XIV) is consistent with experiment. 
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APPENDIX A: DECAY AMPLITUDES OF THREE-BODY B DECAYS 



In this appendix we list the factorizable amplitudes of the 3-body decays B — > 
KKK, KKtt, Ktttt, tttttt. Under the factorization hypothesis, the decay amplitudes are given by 

(PiP2P S \H eS \B) = ^ E aW(p 1 p 2 P 3 |t p |S) j (Al) 

V * p=u,c 

(r) 

where X p = V p bV* r with r = d,s. For KKK and Ktttt modes, r = s and for KKtt and tttttt 
channles, r = d. The Hamiltonian T p has the expression 34] 

T p = aiS pu (ub)v~A ® (su) V -a + a 2 5 pu (sb) V -A ® {uu) V -A + a 3 (sb) V -A ® E(??)v-A 

<? 

+<4^2(qb)v-A <8> {sq) V -A + a 5 (sb) V -A ® E("?<?)v+A 

3 

3 3 
- 2a sE(^) s - p ® 9 e<7 ^ 9 ' )5+p + a 9(^)v-A ® 22 ^e q (qq) V -A 

3 

+ a io z2(Q b W-A ® -e g (sg)y_A, (A2) 

with (<?</) V± A = 97/^(1 =t 75 (qq')s±P = <?(! i 75)9' and a summation over g = u,d,s being 
implied. For the effective Wilson coefficients, we use 

ai ps 0.99 ± 0.037i, a 2 « 0.19 - O.lli, a 3 ps -0.002 + 0.004i, a 5 ps 0.0054 - 0.005i, 
al « -0.03 - 0.02», 4 « -0.04 - 0.008i, a£ ps -0.06 - 0.02*, a c 6 ps -0.06 - 0.006i, 
a 7 ps 0.54 x 10" 4 i, ps (4.5 - 0.5i) x 10 -4 , ps (4.4 - 0.3i) x 10 -4 , (A3) 
a 9 ps -0.010 - 0.0002i, a^ ps (-58.3 + 86.1i) x 10~ 5 , a c w ps (-60.3 + 88.8i) x 10~ 5 , 

for typical at the renormalization scale [i = m^/2 = 2.1 GeV which we are working on. 
Various three-body B decay amplitudes are collected below. 

B -> KKK 



{K°K + K-\T P \B°) = {K + K°\(ub) V -A\B o ){K-\(su) v „ A \0) ai 5 pu + a\ + a p w - (a p + a 



®4 



+(K + K-\(db) v _ A \B°}(K \(sd) V - A \0} (a p 4 - ±a p 10 



+(K \(sb) V -A\B )(K + K-\(uu)v-A\0)(a 2 5 pu + a 3 + a 5 + a 7 + a 9 



+{K°\(sb) v . A \B u )(K + K-\(dd) v ^ A \0) 



1 r \ 
a 3 + a 5 - -(07 + ag) 



+(K"\(sb) v _ A \B")(K + K~\(s S )v^ A \0) 
+(K \sb\B }(K + K-\ss\0){-2a p 6 + a p 8 ) 
+(K+K-\d(l - l5 )b\B°)(K°\s(l + l5 )d\0) (-2a p + a p 8 ) 



a 3 + a p A + a 5 - ~(a 7 + a g + a p l0 ) 
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+(K + K-K \(sd) v _ A \0}(0\(db) v „ A \B }[ a\ - -a% 



+(K + K-K \s l5 d\0}(0\d-f 5 b\B°)(-2al + af 



(A4) 



with rZ = 



2mi 



mi,(/i)(in 2 +rai)(/i) " 



(K + K-K-\T p \B~) = (K + K-\(ub)v-A\B-)(K-\(su) V - A \0} [ ai S pu + a\ + a p w - (a p + a p 8 )r K 

+(K~\{sb) v ^ A \B-){K + K-\{uu) V - A \0)(a 2 5 pu + a 3 + a 5 + a 7 + a 9 ) 



03 + 05 - -jifl-j + ag) 



a 3 + al + a 5 - -(a 7 + a 9 + a p ] 



+(K-\(sb) v _ A \B-){K+K-\(dd)v- A \0) 

+(K-\(sb) v _ A \B-)(K + K-\(ss) v _ A \0) 

+ (K-\sb\B-)(K + K-\ss\0)(-2a p 6 + ag) 

+ (K+^-|n(l - 75 )6|B )(if-|s(l + 75HO) (-2ag + ag) 

+ (K + ^-K-|(su)y_ A |0)(0|(n6)v-A| J B-)f^-ia : 
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+(K+^-K-|s75w|0)(0|w7 5 6|S-)(-2aP + 



(A5) 



Since there are two identical K mesons in this decay, one should take into account the identical 
particle effects. For example, 

{K + K-\{ub) V - A \B-){K-\{su) V - A \Q>) = {K + ( Pl )R- (p 2 )\{ub) v _ A \B~){K~ (p 3 )\(su) v . A \Q) 

+ (K + ( Pl )K-(p 3 )\(ub) V - A \B-)(K-(p 2 )\(su) V - A \0), 

(A6) 

and a factor of \ should be put in the decay rate. 



(K°K K°\T P \B ) = (K°K \(db) v ^ A \B )(K \(sd) v „ A \0)(a p - \a p w - (a p - \a p 

+ (K \(sb) v „ A \W)(K K \(dd) v ^ A \0) 

+ (K \(sb) v „ A \W)(K K \(ss) v „ A \0) 
+(K°\sb\W}(K K°\ss\0}(-2a p + a p ) 

+(K o K o K°\(sd) V - A \0)(0\{db) V - A \B ) ( a p - Ua 7 + a 9 + a 



„K 



03 + «5 - + 09) 



03 + 04 + 05 - 2(07 + ag + af ) 



no, 



+(K K K \s l5 d\0)(0\d l5 b\B°)(-2a p + a p ). 



(A7) 



-rfi 



The second and third terms do not contribute to the purely CP -even decay B — > KsKsKs- 



(K-K"K"\T P \B-) = (K°K \(ub) v ^ A \B~}(K-\(su) v _ A \0) [ ai 8 pu + a p + a p - {a p + a p )r 

^K K~\(db) v „ A \B-)(K°\(- S d) v „ A \0)(a p - \a p w - (ag - \a p )r^) 



+ (K-\(sb) V - A \B-)(K"K"\(dd) v ^ A \0) 



03 + 05 - 2^7 + a 9) 



32 



«3 + <A + a 5 - -(a 7 + a 9 + a p w ) 



+(K-\(sb) V - A \B-}(K K \(ss) v ^ A \0) 

+(K-\sb\B-)(K°K°\ss\0){-2a p 6 + a p ) 
+ (K-K K \(su) v ^ A \0)(0\(ub) v ^ A \B )(a 1 5 pu + a p 4 + a p w 
+(K-K°K \s l5 u\0)(0\u(l - lb )b\B-){2al + 2ag). 
The third and fourth terms do not contribute to the decay B~~ — > K~ KgKg- 

B -> ififTT 



(A8) 



(vr-K+^-|T p |B-) 



(K + K-\{ub) V - A \B-)(Tr-\{du) V - A \0) [ ai 5 pu + a p + a p - {a p + ag)r£ 
+ (7r _ |((i6)v-A|-B _ )(K + i ; r _ |(un)y_. J 4|0}(a2(5p U + a 3 + a 5 + a 7 + a 9 ) 
+(7r-|J6| J B-)( J fC + J fC-|dd|0)(-2a| + ag) 



+(7r-|(J6)y_ A |B-)(K + K-|(s S )y-A|0) 



03 + a 5 - -{a 7 + 09) 



+ (K-\(sb)v-A\B-)(K + 7r-\(d S ) V - A \0)(al - -a 
+ (K-\sb\B-)(K + TT-\ds\0)(-2al + ag) 



+( j PC + K-^-|(Jn)y_ A |0)(0|(u6)y_A|5") (ai<5 pu + + a : 
+(K + j PC-7r-|J 75 n|0)(0|u756|S-)(2a| - ag). 



(A9) 



(K-^+vr-|T p |B-) = (TT + Tr-\(ub) V - A \B-)(K-\{su) V - A \0) L<5 pu + + a p 10 - {a p + ag )r 



J< 



a2$ P u + ~(a 7 + ag) 



+ (^-|(^) y -A|5")(vr + vr-|(t2n)y_ A |0) 
+ (^-|s6|S-)(vr + 7r-|ss|0)(-2ag + ag) 
+(7r-|(J6)y_ A |S-)(^-7r + |(^) y _ A |0)(^ - i< ) 
+ (7r-\db\B-){K-Tr + \sd\0){-2a p (i + ag) 
+ ( J PC-7r + 7r-|(sn)y^|0)(0|(u6)y_ A |S-)(a 1 5 pu + a\ + o? ) 
+ (^-7r + 7r-|s 75 u|0)(0|u75&|5-)(2a^ + 2ag). 



(A10) 



(K°7r + 7r-|T p |5°) 



(7r + 7r-|(J6)y_ A |S°)(if |(^)y_ A |0) 



-a 



10 



(4 



+(A |(s6)v_a|5 )(7r + 7r~|(ttu)v_A|0) 



ct2<5pu + 7: (07 + a 9) 



+ (^|^|5 u )(7r + 7r-|s S |0)(-2a^ + ag) 

+ (7r + |(u6)y_ A |S )(K 7r-|(sn)y_ A |0)(ai + + a?„) 

+ (ir + \ub\B^}(K°iT-\su\0}(-2a p - 2ag) 

+ (Z°7r + 7r-|(^)y_ A |0)(0|(d6)y_A|5 )(ai ( 5 pu + a\ + a p w ) 



+(K°7r + 7r-\s(l + l5 )d\0)(0\d l5 b\B")(2a p 



-f:0> 



(All) 
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(^-vr+vr !^^ ) = {7r + Tr \(ub)v-A\B")(K-\(su) V - A \0) ai 6 pu + al + a 



-Ffl\ 



P 

l 10 



K + 4K 



+ (K-Tr + \(sb)v-A\B°}(ir \(uu)v-A\0) 



+ (ir + \(ub) v ^ A \B )(K-7r°\(s U ) v _ A \0) [ ai 5 pu + a\ + a p w ] 



HA(db) V - A \B Q ){K-Tt+\(- S d) v ^ A \U){al - -aP w ) 

+(ir + \ub\B )(K-7r°\su\0)(-2al - 2af) 
+ (7T°\db\B }(K-TT + \sd\0}(-2al + og) 

+ (^-7r + vr |(^)y^|0)(0|( ( i6)y^|S )K - ±a? ) 
+(^-7r+7r | S -(l + 7 5)d|0)(0|J 7 5fe|^ )(2< - ag). 



(A12) 



(^°7r-7r°|T p |S-) = (7r°7r-|(J6)y_ A |S-)(^ u |(^)y_ A |0) 



no 



(«5 - ^D^x 



+{ICir~\(sb)v- A \B~)(ir \ (uu) V - A | 0) 



+ (7r |(n6)y_A| J B-)(K%-|(su)y_ A |0) [ ai 5 pu + a p A + a p w ] 
+(it-\{db)v- A \B-) (k\° I (sd)v- A 1 0) U 



1 v 
2<o 



+ (vr |u&|£-)(l? ( V|su|0)(-2a£ - 2ag) 

+ {7r-|J6| J B-){i? vr |scZ|0)(-2ag + a|) 

+ (ir 7r~7r°|(^)y_ A |0)(0|(u6)y_ A | J B-)(a 1 «5 pu + a\ + o? ) 

+ (ir°7r-7r°|s(l + 7 5 H0)(0|n 7 5&|5-)(2a£ + 2ag). 



(A13) 



(K\°7r \T p \B ) = (nW\(db) V - A \W)(K"\(sd) V - A \0) 



a2<5p« + 7t(-07 + 09) 



p p 
a 4 ~~ 2 10 



+ (Z 7r |(s6) y _ A |S )(7r |(un)y_ A |0) 
+ (^ |(*)v-A|S )(KV|(^)y„ J 4|0) 

Mn\ \(m) v _ A \0){K°\(sb) v _ A \B°) (a 2 5 pu + 2a 3 + 2a 5 + ^(a 7 + a 9 ) 

+(K°TT \sd\0}(TT \db\B }(-2al + ag) 
+(^V |ss|0)(Z |^|5°)(-2a^ + ag) 
+ (Z 7r vr |(^) y _ A |0)(0|(^) y ^|S )(4 - ±a? ) 



+ (Z U 7r°7r°|5(l + 75 )d|0)(0|(J75fc|5 )(2< - ag). 



(A14) 



— > 7T7r7T 
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(it 7r + 7r \T P \B 



= (ir + ir \{ub) V -A\B ){ir \(du)v-A\fy ai<5 pu + a p A + a p w - (af + af)r£ 

r 3 1 i 

+(-7r~\(db)v-A\B~}(Tr + Tr~\(uu) v ~A\0) \a2$pu -a\ + -(a 7 + a 9 ) + -a^ 

-(7r-|Jfe| J B-)(7r+7r-|Jd|0)(-2og + a|) 

|//|ti;<U|u7 5 6|5~)(2a^ + 2og). (A15) 



+( 



+(tt tt+tt |(du) y _A|0)(0|(n6) V /_ A | J B-)(a 1 5 pu + < 



+ «?o) 



+ (7T~7r + 7T~|d(l + 75 ) 



(7rV+^-|T p |i3 u ) = (^ + 7r°|(ii6)v_ A |B u )<7r-|(dt i )v-A|0) 



(a? + <K 



+(Tr + TT-\(db) v „ A \B°}(iT \(uu) v _ A \0}\a 2 5 pu - a\ + (ag - ^af)r* 



3 

+ -(a 7 + a 9 ) + -a? 



+ ^ + |(u6)y_ A |S u )(7r-7r |(dn)y_ yi |0) [ai5 pu + < + a? ] 

r 3 1 

+(ir°\(db) v ^ A \B ){w + n~\(uu)v-A\fy a 2 5 pu -a p 4 + -(07 + 09) + t^io 



+(7r°|J6|S°)(7r + 7r-|dd|0)(-2a|? + a| 



APPENDIX B: DECAY CONSTANTS, FORM FACTORS AND OTHERS 



(A16) 



In this appendix we collect the numerical values of the decay constants, form factors, CKM 
matrix elements and quark masses needed for the calculations. We first discuss the decay constants 
of the pseudoscalar meson P and the scalar meson S defined by 

{Pip) I <Z27/x75<?i|0) = -ifpPn, {S(p)\q2Jij,qi\0) = fsp», {S\q 2 qi\0) = m s f s , (Bl) 

and {V(p, e)|V^|0) = /ymye* for the vector meson. For the scalar mesons, the vector decay 
constant fs and the scale-dependent scalar decay constant fs are related by equations of motion 

m s 



Vsfs = fs, 



with {is 



(B2) 



miin) - mi(fi) ' 

where rri2 and mi are the running current quark masses. The neutral scalar mesons a, /o and a[j 
cannot be produced via the vector current owing to charge conjugation invariance or conservation 
of vector current: 



fa = ffo = fa° n = 0. 



(B3) 



However, the decay constant fs is non- vanishing. In [48| we have applied the QCD sum rules to 
estimate this quantity. In this work we folow 48f] to use 



4(980) = 460 MeV, f K *(U30) = 550 MeV, 



(B4) 
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at jj, = 2.1 GeV. As for the decay constants of vector mesons, we use (in units of MeV). 



f p = 216, f K . = 218, 



//o(980) = 460, f K * 



550. 



(B5) 



Form factors for B - 
(P(p')\V,\B(p)) 

(Sip'M^Bip)) 
(V{p',e)\V»\B{p)) 
(V(p',e)\A,\B(p)) 



P, S transitions are denned by 41 ] 



(p + p') 



m 



B 



rrii 



in 



B 



nil 



(p + p')n 



in 



B 



m2 



r 



F BS 



(Q 2 ) + 



m, 



nii 



Q„F BS (q 



^uape* v p a p' P V{c 



m B + my 

m B + m v )el l Af v (q 2 

-2my £ -^ qfl [A BV (q 2 



£~ ■ p 

ms + my 



(p + p>)^ V {q 2 ) 



A BV 
^0 



where q = p - p', F x (0) = F (0), ^ 3 (0) = A (0), and 



A 3 (q 2 

where P M = (p + p%, = 
transition in order for the B —. 
quark symmetry 

Various form factors for B 



nip + my 
2niy 



A,{q 2 



nip — my 
r Jjriy 



A 2 (q 2 



(B6) 



(B7) 



(p — p')^. As shown in 65J, a factor of (— i) is needed in B — > S 
-> S form factors to be positive. This also can be checked from heavy 



5* transitions have been evaluated in the relativistic covariant 
light-front quark model 65J]. In this model form factors are first calculated in the spacelike region 
and their momentum dependence is fitted to a 3-parameter form 



F(q 2 



1 — a(q 2 /m 2 B ) + b(q 2 /m 



2 \2 

b) 



(B8) 



The parameters a, b and F(0) are first determined in the spacelike region. This parametrization is 
then analytically continued to the timelike region to determine the physical form factors at q 2 > 0. 
The results relevant for our purposes are summarized in Table IXlXl In practical calculations, we 
shall assign the form factor error to be 0.03. For example, F^FjO) = 0.35 ± 0.03. 

The form factor for B to /o(980) is of order 0.25 at q 2 = [48]. In the qq model for the / (980), 
F B fo = F B f° sin O/y/2. 

For the heavy-flavor independent strong coupling g in HMChPT, we use \g\ = 0.59 db 0.01 ± 0.07 
as extracted from the CLEO measurement of the D* + decay width [39]. The sign is fixed to be 
negative in the quark model 231 ] . 

For the CKM matrix elements, we use the Wolfenstein parameters A 
p = 0.195 and f) 



0.326 



0.806, A = 0.22717, 



521 ] . The corresponding CKM angles are {svo.2(5)ckm 
7 = (59 ± 7)° [521 ] . For the running quark masses we shall use 



0.695 



-0.018 
-0.016 



and 



m b (m b ) = 4.2 GeV, 
m c (mb) = 1.3 GeV, 
m s (2.1GeV) = 90 MeV, 
m d (lGeV) = 6.3 MeV, 



m 6 (2.1GeV) = 4.95 GeV, 
m c (2.1GeV) = 1.51 GeV, 
m s (lGeV) = 119 MeV, 
m u (lGeV) = 3.5 MeV. 



m b (lGeV) = 6.89 GeV, 



(B9) 



The uncertainty of the strange quark mass is specified as m s (2.1 GeV) = 90 ± 20 MeV. 
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TABLE XIX: Form factors of B 
front model 



65] 



7r, K, Kq (1430), p transitions obtained in the covariant light- 



F 


F(0) 




a 


b 


F 


F(0) 




a 


b 


pB-w 


0.25 


1.16 


1.73 


0.95 


F Bn 


0.25 


0.86 


0.84 


0.10 


F BK 


0.35 


2.17 


1.58 


0.68 


F BK 


0.35 


0.80 


0.71 


0.04 




0.26 


0.70 


1.52 


0.64 




0.26 


0.33 


0.44 


0.05 


V b p 


0.27 


0.79 


1.84 


1.28 


A*" 


0.28 


0.76 


1.73 


1.20 


Af p 


0.22 


0.53 


0.95 


0.21 


A*? 


0.20 


0.57 


1.65 


1.05 


yBK* 


0.31 


0.96 


1.79 


1.18 


A BK* 


0.31 


0.87 


1.68 


1.08 


A BK* 


0.26 


0.58 


0.93 


0.19 


A BK* 


0.24 


0.70 


1.63 


0.98 
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